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ABSTRACT 



H_ optimal control theory, based on singular value loop shaping, is used to 
synthesize a controller for the statically unstable longitudinal dynamics of X-29 
aircraft. Two design cases are studied: 2-input 2-output; and 3-input 3-output cases. 
H_ theory provides a direct, effective procedure for synthesizing control laws satisfying 
specified performance objectives and robustness specifications. The 2 I/O case has 
better performance, a faster response and is more robust, than the 3 I/O case. 
Discussion and comparison of results are given. 
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I. INTRODUCTION 



The main task for control system designers is to synthesize a control law which 
maintains system response and error signals to within prespecified tolerances while at 
the same time ensure relatively high robustness. Some measures of system robustness 
commonly used are stability margins, sensitivity and disturbance attenuation. Thus, a 
robust system denotes a system with satisfied system performance, has large stability 
margins, good disturbance rejection and low sensitivity to parameter variations. 

For the single-input single-output (SISO) systems, some techniques, i.e., root- 
locus method, Bode plot and Nyquist diagram, had been well developed to design and 
estimate system performance and robustness. Some of these classical methods can be 
generalized to multi-input multi-output (MIMO) systems, however, not all of them do. 
For instance, from a Nyquist diagram, one can not get a satisfactory notion of 
multivariable stability margins. 

In the past decade, a lot of research on properties of multivariable feedback 
system has been accomplished. Singular value Bode plots of return difference and 
loop gain matrices have emerged as useful indicators of multivariable robustness [Refs. 
1, 2, 3]. Among the techniques that have been recently developed, H_, frequency- 
weighted LQG, and LQG loop transfer recovery optimal synthesis theories have made 
singular value loop shaping a routine matter [Ref. 5]. For a multivariable feedback 
control system, singular value loop shaping involves the manipulation of system loop 
gain over a specified frequency bandwidth to improve system perfonnance and 
robustness. 
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H„ theory provides a direct, reliable procedure for synthesizing a controller. 
While the frequency-weighted LQG optimal synthesis theory (or "H 2 theory") and LQG 
loop transfer recovery theory lead to somewhat less direct but highly effective iterative 
procedures for shaping singular value Bode plots to satisfy singular value loop shaping 
specification [Ref. 5]. 

The aim of this thesis is to present a controller design procedure, bases on the 
H„ theory, for a reduced order, linearized longitudinal dynamics model of the X-29 
aircraft analog backup mode. Some review about the properties of singular value, 
sensitivity reduction and return difference matrices will be given in Chapter II. In 
Chapter 1H, IT. theory and weighting matrix are introduced. Application of H_ theory 
to X-29 longitudinal dynamics model is made in Chapter IV. The design result is 
discussed in Chapter V. 
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D. FEEDBACK PROPERTIES 



For a control system, some properties such as impulse response matrix, transfer 
function matrix and other characterizations of the response to commands can be 
changed by prefiltering the command signal, rather than using feedback loop. But 
some properties, usually termed feedback properties, like stability, sensitivity and 
disturbance attenuation can be altered only through the feedback [Ref. 2]. Here, we 
focus on these feedback properties which also define the robustness of system. 

A. RETURN DIFFERENCE MATRIX 

Consider the standard linear time-invariant multivariable feedback control system 
in Fig. 2.1. It consists of the interconnected plant (G), controller (F) and sensor (H) 
forced by commands (r), measurement noise (n), and plant disturbance (d). The 
prefilter (P) is an optional element used to achieve deliberate command shaping. 

From this configuration, we get 
plant 

y(s) = G(s)u(s) + d(s) (2.1) 

controller 

u(s) = P(s)r(s) - F(s)z(s) (2.2) 

sensor 

z(s) = H(s)(y(s) + n(s>) (2.3) 

where y(s), d(s), n(s), u(s), z(s), r(s) are vectors, and G(s), F(s), H(s) are matrices. 
For the feedback loop, define the transfer matrix 

M(s) = F(s)H(s) 

3 
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Figure 2.1 Standard Feedback Control System 



and matrices 

L 2 (s) = G(s)M(s) 

L,(s) = M(s)G(s) 

are called return ratio matrices [Ref. 2] at the y(s) node and the u(s) node, 
respectively. 

And 

I + L 2 (s) 

I + L,(s) 

are return difference matrices at y(s) node and u(s) node. Similarly, the y(s) node and 
u(s) node inverse return difference matrices are defined, respectively, as 



I + L 2 1 
I + L, 1 



Those terms introduced above are fundamental and important in the control system 
synthesis. Because for a SISO system, the phase and gain margins can be determined 
by the behavior of return difference as a function of frequency. And it turns out, as 
we will see later, the return difference matrix also provides a means of measuring 
robustness in a multivariable system. 

B. SINGULAR VALUES AND MATRIX NORM 

In the design of a feedback loop, one needs to determine the "size" of matrix 
to estimate its properties. For a vector X(x u x 2 , x 3 ,...,xj, the size can be defined as 
the Euclidean norm 

II X II = V x, 2 + x 2 J + x} + + x} 

For a matrix, a more general concept is needed. One way to describe the size of 
matrix A is based on the largest singular value of A, as 

0(A) = max( II AX II / II X II ) 

where II ♦ II denotes Euclidean norms. The singular values a, of a rank n matrix Ae 
C“ are defined as the non-negative square roots of the eigenvalues of A*A, where A* 
is the transpose of the complex conjugate of A. It is convenient to order them as 
follow 

Gj £ a 2 £ <Tj ... £ a„ 

If rank r < n, then there are n - r zeros singular values, that is 

^r+l = ^r+2 = = = 0 
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One uses singular values rather then eigenvalues to describe matrix size since 
eigenvalues can be a very poor indicators of the "size" of matrix [Refs. 5, 14]. The 
singular value is also a good measure of the near-singularity of a matrix, i.e., how 
near the determinant of a matrix comes to being zero. The condition number, the ratio 
of the largest to the smallest singular value, provides information about sensitivity to 
perturbation, which also is an important indicator of robustness. Listed as follows are 
some useful properties of singular value : 

(1) 5(A) = max (II AX II / II X II) 

(2) 5(A) = min (II AX II / II X II) (the least singular value of A) 

(3) 5(A) <, I A,(A) I <, 5(A) , where \ denotes the i-th eigenvalue of A 

(4) If A' 1 exists, 5(A) = 1 / o(A-') 

(5) If A' 1 exists, 5(A) = 1 / 5 (A ') 

(6) 5(AB) £ 5(A) 5(B) 

(7) S° i .,03 = Trace (A‘A) [Ref. 5] 

In SISO system, it is customary to define sensitivity in terms of percentage 

variation in system output resulting from a given percentage change in the plant with 
the input held fixed. A similar way is used to develop the sensitivity function for a 
MEMO system. 

C. SENSITIVITY FUNCTION AND DISTURBANCE ATTENUATION 

Consider a finite dimensional linear time invariant (FDLTI) system in Fig. 2.2, 
where r, d, and n are command, disturbance, and measure noise respectively. If it is 
asymptotically stable, then we have input-output loop error and sensitivity relation as 
follows : 
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Figure 2.2 Finite Dimentional Linear time Invariant System 
input-output 

y = GF(I + GF)-'(r - n) + (I + GF)d (2.4) 

loop error 

e = r - y 

= (I + GF)-’(r - d) + GF(I + GF)’n (2.5) 

system sensitivity 

AHo. = (I + G'Fy'AHoL (2-6) 

where AHc. and AHq,. represent the changes in the close-loop system and nominal 
open loop system caused by changes AG in the plant G. i.e., G' = G + AG [Ref. 1]. 
Here, we define 

S(s) s (I + L(s))' 1 = e/r (2.7) 
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R(s) s F(s)(I + L(s))' 1 = u/r 

T(s) = L(s)(I + L(s))' 1 = I - S(s) = y/r 



( 2 . 8 ) 



(2.9) 

where L(s) = G(s)F(s) is the loop transfer function, S(s) is known as the sensitivity 
function matrix, T(s) is complementary sensitivity function. It can be seen that loop 
transfer function L(s) determines the other three matrices. And, as will be seen, the 
singular value bode plots of these matrices play important role in the design of 
multivariable robust control system. 

From Eq. (2.5), it is seen that the loop error " e " resulting from disturbance 
inputs can be made small by making the sensitivity matrix, or inverse of the return 
difference matrix (I + GF)' 1 , small. From Eq. (2.6) the system sensitivity to plant 
variations may be also reduced by decreasing the sensitivity function. One can 
decrease the size of the sensitivity function by making 

G(S(j(D)) = o[(I + L(jco))-'] 
small, or (from singular value properties (5)) 

5(1 + L(j(0)) (2.10) 

large. 

Note also 



G(I + L(ja>)) = 5(L(jco)) for 5(L(jo))) » 1 
This tells us that large loop gains or " tight " loops yield good performance. But one 
can not make loop gains arbitrarily high over arbitrarily large frequency. Since we can 
see, from Eq. (2.5), large 5(L(jco)) over a large frequency range not only make errors 
due to r and d small, but make error due to n large. Large gains can also lead to 
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instability. There are certain performance tradeoffs between command / disturbance 
error reduction and sensor noise error reduction and stability to be satisfied. 

D. STABILITY AND STABILITY MARGIN 

In SISO system, gain margin and phase margin are the stability margins that are 
commonly used. Both are conveniently expressed in terms of the magnitude of the 
return difference, I + L{j(o). In MIMO system, the return difference retains its 
importance, but the concept of gain margin and phase margin become problematic. 

The stability of the system in Fig. 2.2 is directly related to the matrices S(s), 
R(s) and T(s). To find the stability margin for the system is to determine a lower 
bound on the " size " of smallest perturbation to the plant that will destabilize the 
system. And the singular value Bode plot of S(s) and T(s) can be used as the measure 
of stability margins. 
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III. H_ CONTROL DESIGN 



IT. control theory provides a direct, reliable procedure for synthesizing a feedback 
controller designed to meet singular value loop shaping requirements. Its capability 
includes addressing the full range of stability margins, sensitivity, and robustness 
optimization and response-shaping problems that can be formulated within the singular 
value Bode plot framework. The standard configuration of an H„ problem is shown 
in Fig. 3.1. The design objective is to fmd a stabilizing controller F(s) for the 
augmented plant P(s) while keeping the robustness specifications satisfied. 




Figure 3.1 Standard H„ Small Gain Problem 
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The purpose of this chapter is to present a brief overview of H„ optimal theory 
and a control system design procedure based on this theory. 

A. CONTROL SYSTEM FORMULATION 

As shown by Eq. (2.5), (2.6), and (2.10), the loop error due to disturbance and 
system sensitivity to plant variations can be reduced by suppressing the sensitivity 
function S(s). (There is also some stability requirements desired in the robust system 
synthesizing.) These goals can be achieved by properly selecting the weighting 
functions in the design procedure. 

Shown in Fig. 3.2, is the compensation configuration of system in Fig. 2.2, where 
P(s) is the augmented plant with weighting functions W t , W 2 and W 3 which penalizing 
the error signal, control signal and output signal respectively. 



Augmented Plant P(s) 




Figure 3.2 Compensated System with Augmented Plant P(s) 
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From this Figure, we get 
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Where the augmented plant with weighting functions is denoted as 
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A state space realization of P(s) is given by 
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( 3 . 2 ) 
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If the weighting transfer function is improper (i.e., has more zeros than poles), as can 
be the case for W 3 , then no state space realization exists. The produce W 3 G is in 
this case still proper and can be represented by a state space realization. Some 
calculations are need to compute the state space realization of the augmented plant. 
From Fig. 3.1 one can define the transfer function T ylul as 

y, = . 

This transfer function can be expressed in terms of the weighting matrix and previous 
defined loop matrices, i.e., 



'yu ' 




"w,s " 


y.b 


= 


W : R 


_y.c . 




W 3 T 



And the close loop transfer function from input u t to output y t is denoted as 



T 



ylul 



w,s 

W : R 

W 3 T 



Where S, R, and T are loop matrices defined in Eqs. (2.7), (2. 8), (2.9). 

The compensated system is said to be internally stable if the augmented A matrix 
of the compensated system is stable, i.e., when external input signal u, equals to zero, 
the states of both P(s) and F(s) will go asymptotically to zero for any initial 
conditions. And the controller is said to be stabilizing. 
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B. H„ OPTIMAL CONTROL PROBLEM 

The H„ optimal control problem is defined as the following small gain 
problem: 

For the given transfer function matrix P(s), find a stabilizing controller F(s) 
such that the close loop transfer function matrix T yIuI is internally stable and its 
infinity-norm is less than or equal to 1, i.e., 

II T ylul II < 1. (3.9) 

Where the H„ - norm of matrix A is defined as 

II A IL = sup a(A(jco)). (sup : the least upper bound) 

This problem is called a " small gain problem " because the closed loop gain is small 
(i.e., less than or equal to one). As mentioned before, the largest singular values of 
S(s) determine the plant disturbance attenuation, since S(s) is the same as the close 
loop transfer function from disturbance d to plant output y (Fig. 2.2). One may 
specify the system disturbance attenuation performance with the frequency-dependent 
weighting function W„ that is 

a(S(jco)) £ I Wf'doo) I 

Now, consider the system in Fig. 3.4, where A a (s) and A M represent the additive 
plant perturbations and the multiplicative plant perturbations respectively. Defining the 
size " of A A (jco) and A M (jco) as a(A A (jw)) and A M < j<o)). the stability theorems are 
obtained as follows: 
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PERTURBED PLANT 




Figure 3.3 Additive / Multiplicative uncertainty 
Robustness Theorem 1: 

If the system in Fig. 3.3 is stable as both A a and A M are zero. Let A a = 
0, the size of the smallest stable A Js) for which the system become unstable 

is 

1 

o(A(jo))) = 

o(T(jo))) 

Where T(jco) is the complementary sensitivity function of the system. It is seen that 
decreasing the value of d(T(j(o)) will increase the size of A(j(o), that is, increase the 
stability margin. 

Robustness Theorem 2 : 
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Suppose the system in Fig. 3.3 is stable when A a and A M are both zero. Let 
A m = 0. Then the size of the smallest stable A A (s) which destabilizes the 
system is 

1 

a(A A (jco)) = — 

a(R(jco)) [Ref. 5] 

We can see that the smaller d(R(jco)) is, the greater will be the stability margin. 

With properly selected weighting function matrices W 2 (jco) and W 3 (ja>), the 

stability margins of the control system can be specified by the following singular value 

inequalities, 

d(R(jco)) < I W 2 '(jco) I 
a(T(j<o)) < I Wj-'Cjto) I 

Together with the performance specification W,(jco), the design specifications of the 
robust control system may be written as 



o(S(jco)) < 1 


W,-’(j®) 1 


(3.3) 


a(R(jco)) < 


1 W 2 '(jco) 1 


(3.4) 


VI 

3 

‘C? 

H 

ID 


1 W/'(jco) 1. 


(3.5) 



Thus, for a robust control system synthesis, the requirements may be presented as 
reducing the plant disturbance effect as much as possible while meeting the control 
inputs and the robustness constraints. This goal can be achieved by the singular value 
loop shaping iterative procedure that will be introduced in the following sections. 

As a consequence of singular value properties and the definition of H„ norm, 
these singular value inequalities can be combined into the single infinity norm 
specification in terms of the close loop transfer function T ylul of the form 
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where 



II T yUI IL ^ 1 



T 



ylul 



W,S 

w 2 r 

W,T 



Note, the matrix T ylul is partitioned into submatrices representing performance 
requirements and stability constraints. 



C. WEIGHTING MATRIX SELECTION 

Since the plant disturbance can be attenuated by reducing the value of d(s(jco)) 
continuously, a constant scaler weighting gain y is added to W, to facilitate the 
singular value looping shaping process. That is 

a(S(jco)) £ I Y 1 * W^co) I 

So, as we increase the value of y, the sensitivity function is suppressed down 
continuously. On the other hand, from Eq. (2.9) (T(s) + S(s) = I), we know that T(s) 
goes to I as y increases (since S(s) goes to zero). This means the complementary 
sensitivity is forced against its upper constraint W 3 ''(s) as the value of y increases. 
The closed loop bandwidth, therefore, is widened. This improves the system 
performance, i.e., make system response faster. At the same time, the singular values 
of closed loop transfer function T ylul is pushed up to upper limit. The singular value 
loop shaping process stops as the condition in eq. (3.9) is reached. 

It is interesting to note that, 

S(s) = (I + L(s))' 1 = L(s)' 1 
as 2(L(jco)) » 1, ( i.e., the system has large loop gains) 
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1 / a(S(jw)) « s(L(j(o» 



i.e., 

Similarly, as a(L(j(o)) « 1, 

T(s) = L(I + L(s))' 1 « L(s) 

that is, 

a(T(j©)) « a(L(jco)) 

These can be seen in Fig. 3.4, i.e., above the zero dB line, S(L(j(o)) » 1, 1 / a(S(jot>)) 
= s(L(jco)), while below the zero dB line, a(L(ja>)) « 1, a(T(j(o)) = a( L(ja>)). Thus 
the robust system design consists of the use of high loop gains at low frequencies 
where the performance specifications are defined and the use of low loop gains at high 
frequencies where the robustness constraints lie [Refs. 5, 13]. The system performance 
will be limited within the range between a(G(s)F(s)) curve and g(G(s)F(s)) curve. 

The weightings W, and W 3 are selected as transfer function matrices of which 
diagonal elements are frequency-dependent constants. And the size of these weighting 
matrices is the same as the number of plant output states. This small gain problem 
is solved by using a software program, i.e., " hinf.m " provided by MATLAB 
ROBUST-CONTROL TOOLBOX package. The number of states of the stabilizing 
controller F(s) produced by " hinf.m " is the same as that of P(s). In addition it is 
required by " hinf " that the D 12 submatrix of P(s) (see Eq. (3.1)) is full column rank. 
An easy way to ensure that this is the case is to choose the weighting matrix W 2 with 
a invertible " D-matrix ", e.g., W 2 (s) = el, where £ is any non-zero number [Ref. 5]. 
The parameters considered in selecting these weightings are gain, order, attenuation 
and comer frequency, which are used to achieve small a(S(s)), i.e., large loop gains, 
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Figure 3.4 Singular Value Spec’s on S(s) and T(s) 

over the broadest possible bandwidth subject to the robustness constraints. And it 
should be noted that the controller F(s) depends on the weighting matrices selected. 
Once the weightings are selected, the control system design is in a computer-aided 
design environment and requires only " one parameter " y-iteration. The flow chart 
of the H_ iterative procedure is shown in Fig. 3.5. 
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Figure 3.5 H„ Singular Value Loop Shaping Iterative Process 
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IV. H_ CONTROLLER SYNTHESIS FOR X-29 FIGHTER 

The purpose of this Chapter is to present a modem aircraft controller design 
procedure based on the H_ optimal control theory. The model used is the longitudinal 
dynamic model of X-29 fighter. The X-29 is a single seat forward swept wing (FSW) 
demonstrator aircraft build by Grumman Co., which is designed to be a new generation 
of tactical fighter with the advantages of light weight, low cost and high efficiency. 
The computer-aided design program, MATLAB-ROBUST CONTROL TOOL BOX, is 
used in controller synthesis. The brief plant description and design objective are given 
in the following sections. The design procedure and results are also presented. The 
script files specifically written or modified for this problem are listed in appendix A. 

A. PLANT DESCRIPTION 

The aerodynamic advantages of forward swept wing design include improved 
maneuverability, with virtually spin-proof characteristics, better low-speed handling and 
reduced stalling speed. The FSW X-29 (see Fig. 4.1) was selected primarily because 
its multiple, independently controlled surface make it an ideal candidate for multiloop 
synthesis of advanced control mode. 

The aircraft is designed to have a -35% stability margin at subsonic speed. A 
wing strake extends aft from the trailing-edge at each wing root, each strake with a 
trailing-edge flap which has its own integrated servo actuator to augment foreplanes 
for pitch control. All-moving canard surfaces, one on each side of the center fuselage, 
outboard of engine inlet ducts and operated by servo actuators for primary pitch 
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Flap 





Figure 4.1 FSW X-29 Fighter 

control. In addition the aircraft is equipped with flaperons [Ref. 10]. 

The X-29 longitudinal dynamic model considered in this study is an analog 
reversion mode with the aircraft trimmed at 0.5 mark, 30,000 feet. The design point 
was one of the critical nodes in the analysis of the aircraft. The 14-state reduced- 
order model includes a short period approximation of the aircraft longitudinal dynamics, 
vertical velocity w and pitch rate q, and 3 fourth order actuator dynamics for the 
longitudinal control surfaces respectively, i.e., the canards, flaperons and strakes. 

For the first design, a 2-input, 2-output configuration is used. Two separated 
commands r, and r 2 are input to the three control surface actuators with r, controlling 
the canards and strakes and r 2 controlling the flaps. The outputs of the system are two 
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aircraft states, w and q. The physical configuration of this 2-input, 2-output open loop 
actuator/aircraft dynamics model is shown in Fig. 4.2. 

The second design example considered is a 3-input, 3-output controller synthesis. 
Three separated commands r,, r 2 and r 3 are input to the three control surfaces canards, 
flaps and strakes respectively. The output to be controlled are w, q and the canard 
control input 5 C . The physical configuration for this model is shown in Fig. 4.3. 
Comparison is made between this two proceeding designs. 

The state variables are listed in Table 4.1. The state space realization of the 
open loop plant model G is listed in Appendix B. The poles of this system are listed 
in Table 4.2. It is seen that there is a positive pole, i.e., 1.9550, on the real axis. 
That means the X-29 has an unstable short period mode. 

The 14-state model was scaled to improve the numerical conditioning of its state 
space representation. The w state was transformed to angle of attack a using the initial 
forward velocity u^, i.e., a = w / u 0 . The units of the actuator third derivative states 
were transformed from rad/sec 3 to le+04 rad/sec 3 . This scaling was effective in 
reducing the condition number of the system state matrix A g from an order of 
magnitude of 10 10 to 10 4 [Ref. 13]. 

The singular value plot of the uncompensated plant G(s) is presented in Fig. 4.4. 
The solid curve and dashed curve represent the maximum and minimum singular values 
of G(s) respectively. As can be seen, at low frequency, the open loop plant has poor 
disturbance attenuation, high sensitivity to plant variations and modeling errors, and a 
small control bandwidth. These properties will be improved by singular value loop 
shaping procedure presented in the following sections. 
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Figure 4.2 X-29 2-Input, 2-Output Open Loop Configuration 
Source : W. L. Rogers [Ref. 13] 
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Figure 4.3 X-29 3 I/O Open Loop Configuration 

Source : W. L. Rogers [Ref. 13] 
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Table 4.1 Discretion of The Uncompensated X-29 state variables 



State 


Description 


Units 


a 


angle-of -attack 


rad 


q 


pitch rate 


rad/sec 


5 c 


canard control input 


rad 




flap control input 


rad 




strake control input 


rad 




canard control rate 


rad/sec 




flap control rate 


rad/sec 




strake control rate 


rad/sec 




canard control accel . 


p 

rad /sec 


h 


flap control accel. 


rad/sec 2 




strake control accel. 


rad/sec 2 


*e 


canard control jerk 


le+04 rad/sec 




flap control jerk 


le+04 rad/sec 




strake control jerk 


le+04 rad/sec' 



Source : W. L. Rogers [Ref. 13] 
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Table 4.2 Uncompensated X- 29 open loop poles 



poleg - 



2 . 27 4 6e+02 + 
2.274 6e+02 - 
1 . 4455e+02 
1 . 4491e+02 
1. 9550e+00 
1 . 003 le+02 
2 . 7 155e+00 
5 . 2518e+01 + 
5 . 2518e+01 - 
5.2506e+01 + 
5 . 2506e+01 - 
5. 0067e+01 
2 . 0172e+01 
2 . 0115e+01 



2 . 3201e+02i 
2 . 3201e+02i 



4 . 8255e+0li 
4 . 8255e+0li 
4 . 8410e+0li 
4.8410e+0li 



X-29 SV PLOT OF PLANT G 




FREQUENCY - r*d/sec 



Figure 4.4 Singular Value Plot of The Uncompensated Plant G(s) 
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B. DESIGN SPECIFICATIONS AND WEIGHTING FUNCTIONS 



Since the design objective is to attenuate the disturbance effect as much as 
possible within loop gain constraint and robustness specifications. The design 
specifications, in this case, are : 

Performance Spec.: 

Minimize the sensitivity function as much as possible. 

Robustness Spec.: 



1) . Attenuate the close-loop singular values of the complementary 

sensitivity matrix by 20 db at frequencies beyond a> = 100 rad/sec 
(15.9 Hz) (This will ensure that the system will have sufficient stability 
margin to tolerate variations the loop transfer function of magnitude as 
large as a factor of ten at frequency (0 = 100 rad/sec). 

2) . a second-order roll-off beyond (0 = 100 rad/sec (This ensures the 

controller is proper and matches the uncontrolled drop-off). 



With these considerations, we have covered the most important issue in feedback 
control system design. Once the design specifications are decided, the selection of 
weighting matrices is the next step. For the 2-input 2-output system, the weighting 
matrices are chose as follows : 



.01 ( 1 + s / .01 ) 



(y w,(s))-' = r 1 * 



* I 



(4.1) 



( 1 + s / 100 ) 



(2 x 2) 



W 2 (s) = - 0.025 * I. 



(4.2) 



(4 x 4) 
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W,-'(s) = 



1000 



(4.3) 



* I 

s 1 (2 x 2) 

to penalize the error " e ", control input " u " and output " y " respectively. 

The singular value plot of W,'‘(s), W 3 ''(s) weighting functions is given in Fig. 
4.5. The robustness specification W 3 ''(s) has 20 dB drop-off and 40 dB roll-off at 
0) = 100 rad/sec which satisfied the design requirements. As mentioned, W,''(s) and 
W 3 ''(s) should be so chose to ensure the 0 dB crossover frequency of W, is sufficiently 
lower than the 0 dB crossover frequency of W 3 ''(s). W 2 (s) is selected to penalize the 



X-29 Design Specifications 




Figure 4.5 Singular Value Plot of Weighting Functions 
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control input u, and the weighting factor 0.025 is chosen to confine the actuators’ 
response within their design limits. The selected W 2 (s) also ensures that the submatrix 
D [2 of P(s) has full column rank as required in the two Riccati equation solution to the 
H„ problem. The FL theory requires that D 12 (s) matrix is full column rank in order 
for an optimal H„ controller to eliminate infinite energy inputs. If we decrease the 
weighting factor, more energy will input to the system, in this case, the actuators’ 
response may exceed the limits [Ref. 13]. 

C. DESIGN PROCEDURE 

Before starting the singular value loop shaping procedure, it is necessary to find 
the balanced state space realization for the unaugmented plant for a better numerical 
condition. This can be seen from the different condition numbers of state matrix A g 
before and after balanced realization, which are 5.526e+04 and 2.5567e+02 respectively. 

1. 2-INPUT, 2-OUTPUT DESIGN CASE 

With the selected weighting matrices, the augmented plant P(s) of the 
balanced plant is of order eighteenth, W,(s) and W 2 each adding two states to the 
system G(s). Since W 3 (s) is not a proper transfer function Eq. (4.3) (i.e., it has more 
zeros than poles), so it has no state space realization. But W 3 (s)G(s) is proper and can 
be realized in state space form. 

With the computer-aided design program, y is chosen as 1 for the first try 
of the iteration process. As continuously increasing the value of y, the singular value 
of cost function T yUI is pushed up to its upper limit gradually until no solution exists 
for any larger y. After several iterations, y equals to 2.6, and the maximum singular 
values of T y , u1 reached the all pass limit (i.e., 0 dB line). Fig. 4.6 is the singular value 
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COST FUNCTION Tylul (Gamma = 2.6) 




Frequency - rad/sec 

Figure 4.6 Cost Function T y)o! for y = 2.6 

plot of cost function T ylo) , the solid line represents the maximum singular values and 
dashed line represents the minimum singular values. Fig. 4.7 and 4.8 are the singular 
value plots of complementary sensitivity and sensitivity, the dashed line and dotted line 
again represent the maximum and minimum singular values respectively. The 
complementary function is pushed flat against the specified limit W,' 1 , as the value 
of gamma is increased, the closed loop bandwidth is widened (to 30 rad/sec) and the 
response speed is increased. Similarly, the sensitivity function is pushed down against 
its limit y'W ,' 1 as y value increased. The singular value of I + G(s)F(s) loop gains 
are increased to about 50 dB at low frequency for this design, which markedly 
improves the X-29 performance properties. This large loop gains indicate good 
disturbance attenuation and low sensitivity to plant variations over a control bandwidth 
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COMP. SENSITIVITY FUNCTION AND 1/W3 




Frequency - rad/sec 



Figure 4.7 Complementary Sensitivity Function T(s) and W 3 '' for y = 2.6 

SENSITIVITY FUNCTION AND l/Wl 




Frequency - r*d/»ec 

Figure 4.8 Sensitivity Function S(s) and W,' 1 for y = 2.6 
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of about 3 rad/sec and satisfy the stability requirements. Note, in Fig. 4.8, W,''(s) = 
a(S(jco)) at low frequency. And there is about 40 dB drop-off from to = 10 rad/sec to 
co = 0.1 rad/sec which is the range of interesting. The state space realization of the 
18-state stabilizing controller F(s) is listed in Appendix C. Two uncontrollable states 
are removed using PRO-MATLAB function " minreal.m ". A fast transient poles is 
eliminated while retaining the most important characteristics of the closed loop system. 
Thus, the final controller has 15 states. The poles of this reduced-order controller are 
listed in Table 4.3. 



Table 4.3 Poles of The Reduced-Order Controller 



poleacpr = 

-2 . 1954e+02 

-6 . 5170e+01 + 9 . 5444e+01i 
-6. 5170e+01 - 9 . 5444e+01i 
-1 . 4487e+02 
-1. 3577e+02 

-5 . 3024e+01 + 4.9271e+0li 
-5. 3024e+01 - 4.9271e+01i 
-4 . 1572e+01 + 3.1587e+01i 
-4 . 1572e+01 - 3.1587e+01i 
-1 . 3180e+01 + 1. 9602e+01i 
-1 . 3180e+01 - 1 . 9602e+01i 
-3 . 0118e+01 + 1 . 8361e+01i 
-3 . 0118e+01 - 1 . 8361e+01i 
-1. 0000e-02 
-1. 0000e-02 
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Figure 4.9 X-29 Feedback Configuration 

Fig. 4.9 shows the closed loop system of compensated X-29, where F(s) is 
the 15th order controller and G(s) is the 14th order X-29 plant. Since F(s) is placed 
in series with plant G(s), the command vector r is composed of the reference 
commands r, and r 2 to controlled outputs a and q. So , the closed loop compensated 
system has 2 inputs, 2 outputs and 29 states. The poles of this reduced-order closed 
loop system are listed in Table 4.4. The singular value plot of output, inverse-return 
difference matrix I + (G(s)F(s))' 1 (also called the output, multiplicative return difference 
matrix) for the compensated system is shown in Fig. 4.10. The minimum singular 
value of this matrix, i.e., a rain [I + (G(ja>)F(jco))'’l. provides a measure of the aircraft's 
gain and phase margins with respect to multiplicative modeling errors using the 
universal gain and phase margin curve (Fig.4.11) [Ref. 4], 
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Table 4.4 Poles of The 2 I/O Closed Loop System 



-2 . 2739e+02 
-2 . 2739e+02 
■2 . 1723e+02 
-6 . 487 6e+01 
■6 . 4876e+01 
■1. 3558e+02 
■1. 4595e+02 
■1. 4469e+02 
■1. 4469e+02 
■1. 1118e+02 
■4 . 9782e+01 
■4 . 9782e+01 
■5. 2188e+01 
•5. 2188e + 01 
■5 . 24 64e+01 
■5. 2464e+01 
■4. 1968e+01 
•4 . 1968e+01 
•4 . 1045e+01 
2 . 7932e+01 
2 . 7932e+01 
1. 1994e+01 
1. 1994e+01 
1 . 9603e+00 
2 . 6697e+00 
3 . 8935e+00 
1.0123e+01 
1 . 6232e+01 
1 . 9727e+01 



+ 2 . 3194e+02i 

- 2 . 3194e+02i 

+ 9 . 6310e+01i 

- 9 . 6310e+01i 



+ 9 . 2844e-02i 

- 9 . 2844e-02i 

+ 4 . 9587e+0li 

- 4 . 9587e+01i 
+ 4 . 8038e+0li 

- 4 . 8 03 8e+0li 
+ 4 . 8363e+01i 

- 4.8363e+0li 
+ 3 . 2086e+0li 

- 3 . 2 08 6e+01i 

+ 1 . 3447e+01i 

- 1 . 3447e+0li 
+ 1 . 1959e+01i 

- 1 . 1959e+01i 



It is seen in Fig. 4.10 that + (G(s)F(s))' 1 ] drops to about -3 dB at 

frequencies between 0.3 rad/sec and 20 rad/sec. Entering the vertical axis of the 
universal curve at singular value of 0.71 (i.e., -3 dB). The gain and phase margins of 
compensated X-29 near 0 dB crossover frequency are -1 1 dB to +5 dB and ±42 deg 
respectively. This is more stable than the -8 dB to +4 dB and ±35 deg gain and phase 
margins for the typical design of a fighter aircraft. 
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,„(I+G(ju)) a.i„(I+Q' , ( ju») ) 



X-29 SV PLOT (I + inv(GF)) 




FREQUENCY - rad/sec 

Figure 4.10 Singular Value Plot of The Output Inverse-Return Difference 
Matrix for 2 I/O System 




Figure 4.11 Universal Gain and Phase Margin Curve 
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Shown in Fig. 4.12 are the closed loop Bode plots of a and q respectively. 
The solid curves are the response to r, and dashed curves are to r 2 . It is seen that, at 
low frequency, a is "all pass" (i.e., has 0 dB gain) to input 1 but has very small 
response (-50 dB) to input 2. Similarly, as seen that q is "all pass" to r 2 but has small 
response to r,. This indicates the decoupling of a from q to r, and r 2 . Note, the 
system performance at high frequencies satisfy the robustness requirements for both 
cases. 

The time response of a and q to 1 deg/sec (0.01745 rad/sec) step input are 
presented in Fig. 4.13 and 4.14. It is seen from Fig. 4.13 that the response of a is fast 
while with rise time of 0.3 seconds and 0.1 second to rj and r 2 respectively. We also 
can see that the decoupling of a and q to the two commands from both plots. If we 
reduce the weighting function W 2 (s), more energy will be input to system, and that 
energy increases the response speed and the decoupling will be more pronounced. 

Referring to Fig. 4.15 and 4.16, the time response of control deflections 5 C , 
5„ and 5,. It is seen that the peak magnitudes are 0.013 rad, 0.006 rad, and 0.007 rad 
for 5 e , 5 f , and 6, respectively, which are well within the X-29’s control surface 
deflection limits listed below : 

(1) . canards (leading edge) : 30 deg up / 60 deg down, 

(2) . flaps (trailing edge) : 10 deg up / 25 deg down, 

(3) . strakes (trailing edge) : 30 deg up and down. 

The peak control rates (see Fig. 4.17 and 4.18) are 0.0013 rad/sec, 0.0004 rad/sec, and 
0.0008 rad/sec respectively, which are much smaller than the X-29 actuator minimum 
design requirements of 
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NIVO ap - NIVO 



X— 29 CLOSED LOOP BODE PLOT OF " ALPHA " 




FREQUENCY - rad/sec 



X-29 CLOSED LOOP BODE PLOT OF " q " 




FREQUENCY - rad/sec 

Figure 4.12 Closed Loop Bode Plot of a and q for 2 I/O System 
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TIME - sec 



Figure 4.13 Step Response of a for 2 I/O System 



RESPONSE OF " q ” TO 0.01745 rad / 1 sec STEP INPUT 




TIME - sec 



Figure 4.14 Step Response of q for 2 I/O System 
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X-29 DC FOR 0.01745 rad / 1 sec STEP INPUT 




TIME - sec 



xlO- 3 



X-29 DF FOR 0.01745 rad / 1 sec STEP INPUT 
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Figure 4.15 Step Response of Canard Control Input 5 C and Flap Control 
Input 5 f 
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110-3 X— 29 DS FOR 0.01745 rad / 1 sec STEP INPUT 




Figure 4.16 Step Response of Strake Control Input 5, 



xio- 3 X-29 DCDOT FOR 0.01745 rad / 1 sec STEP INPUT 




TIME - aec 

Figure 4.17 Step Response of Canard Control Rate S c 
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STRAKE CONTROL RATE - rad/sec FLAP CONTROL RATE - rad/sec 



xlO -4 X-29 DFDOT FOR 0.01745 rad / 1 sec STEP INPUT 




xlO -4 X-29 DSDOT FOR 0.01745 rad / 1 sec STEP INPUT 




Figure 4.18 Step Response of 5 f and 5, 
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(1) . canards : 1.75 rad/sec 

(2) . flaps : 0.87 rad/sec 

(3) . strakes : 0.52 rad/sec. 

The response is well within the actuator limits but one also has a very small step- 
input. 

2. 3-INPUT, 3-OUTPUT DESIGN CASE 

The second design example we use is the same X-29 dynamic model as in 
the previous example, but with different input and output vectors, i.e., r„ r 2 , and r 3 to 
control canards, flaps and strakes separately and outputs are a, q, and 5 C (see Fig. 
4.3). Thus, it is a 3-input, 3-output system with 14 states. 

The same design specifications are used for this case, and the weighting 
functions are chose as follows : 

.01 ( 1 + s / .01 ) 

(YW,(s))-' = y 1 * * I (4.4) 

( 1 + s / 100 ) (3x3) 

W 2 (s) = -0.025 * I. (4.5) 

(6 x 6) 

1000 

W/'(s) = * I (4.6) 

s 2 (3 x 3) 

The augmented plant P(s) is an 20th order system with W,(s) and W 2 (s) each adding 
three states to the plant G(s). Following the same iterative procedure as of the 
previous example, the cost function reaches its all pass limit as y = 1.3 (see Fig. 4.19), 
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COST FUNCTION Tylul (Gamma = 1.3) 




then the iterative process stops. The solid, dashed, and dotted curves in this plot 
represent the maximum, medium, and minimum singular values of T ylu , respectively. 
The controller F(s) produced has 20 states, that is the same as number of states of 
P(s). Three states are removed using " minreal.m " function and one state is 
eliminated by " obaLreal.m " program. So, the final controller has 16 states, the 
controller / plant series has 30 states in total (uncompensated plant G(s) has 14 states). 
The poles of this reduced-order closed loop system are listed in Appendix E. 
Comparing Fig. 4.20 with Fig. 4.7 and 4.8 shows that this 3 I/O system are not as 
robust as the 2 I/O model. The 3 I/O system has smaller disturbance attenuation, 
larger sensitivity to plant variations and modeling errors, smaller control bandwidth and 
closed loop bandwidth. Fig. 4.21 is the singular value plot of output, inverse-return 
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SINGULAR VALUE 



COUP. SENSITIVITY FUNCTION AND 1/W3 




Frequency - rad/sec 
SENSITIVITY FUNCTION AND 1/Wl 




Frequency - rad/sec 

Figure 4.20 Complementary Sensitivity T(s) and Sensitivity S(s) for y = 1.3 
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X-29 SV PLOT (1 + inv(GF)) 




FREQUENCY - rad/sec 

Figure 4.21 Output, Inverse Return Difference Matrix for 3 I/O System 

difference matrix. It is seen that the a „ J I + ((G(s)F(s)) 1 ] (i.e., the dotted curve) 
drops to approximately -6 dB between 0.2 rad/sec and 13 rad/sec. Enter the vertical 
axis of the universal curve (Fig. 4.11) at singular value of 0.5 (i.e., -6 dB), the gain 
and phase margins near 0 dB crossover frequency are -6 dB to +3 db and ± 28 deg 
respectively. It indicates that the stability margins of the 3 I/O system, with respect 
to multiplicative modeling errors, are less than the values desired of a fighter aircraft, 
i.e., -8 dB to +4 dB, ±35 deg gain and phase margin respectively. Obviously, it is 
also not as stable as the 2 I/O system discussed before. 

From Fig. 4.22, the Bode plots of a and q to the 3 inputs, we also can see 
that oc is decoupled from q to r ( and r 2 . And it is interesting to see that both cc and 
q have very little response to r,. This also can be seen from Fig. 4.23, i.e., r 3 has 
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GAIN - dB GAIN 



X-29 CLOSED LOOP BODE PLOT OF " alpha 




FREQUENCY - rad/sec 



X-29 CLOSED LOOP BODE PLOT OF ” q " 




FREQUENCY - rad/sec 

Figure 4.22 Closed Loop Bode Plots of a and q for 3 I/O System 
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Response of Alpha to 0.01745 rad/sec Step Input 




Time - sec 




0 0.5 1 1.5 2 2.5 3 3.5 4 



Time - »cc 

Figure 4.23 Step Response of a and q for 3 I/O System 
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almost no effect on both a and q. This indicates that there is not much difference 
to system response between 2 inputs design or 3 inputs design. 

Shown in Fig. 4.24 and 4.25, the peak values of control deflections 5 C , 5 f , 
and 5, (i.e., 0.005 rad, 0.0017 rad, and 0.015 rad respectively) are all well within the 
design limits listed before, and are also smaller than those of the 2 I/O system since 
the loop gains are smaller. Similar results are obtained from the plots of control rates 
(Fig. 4.26 and 4.27), i.e., all responses stay within their design limits. 

Thus, we may conclude that the performance of this 3 I/O control system satisfies 
the design specifications but is not better than those of 2 I/O system designed 
previously. 
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Flap Input - rad Canard Input - rad 



x!0‘ 3 



Response of DC to 0.01745 rad/sec Step Input 




TIME - sec 



xl0~ 3 Response of DF to 0.01745 rad/sec Step Input 




Time - »ec 

Figure 4.24 Step Response of 5 C and 5 f for 3 I/O System 
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Response of DS to 0,01745 rad/sec Step Input 




Time - sec 

Figure 4.25 Step Response of 8, for 3 I/O System 

xl0~ 4 Response of DCDOT to 0.01745 rad/sec Step Input 




Time - sec 

Figure 4.26 Step Response of 5,. for 3 I/O System 
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Strake Control Rate - rad/sec Flap Contro , Rate _ rad/ , ec 



xl0~ 3 Response of DFDOT to 0.01745 rad/sec Step Input 




Time - sec 



rlO -3 Response of DSDOT to 0.01745 rad/sec Step input 




Time - sec 

Figure 4.27 Step Response of S, and 6, for 3 I/O System 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Although the mathematic work behind the H„ theory is very complicated, the 
implementation of this method is easy using a well developed computer-aided design 
software. From the design procedure presented above, it is seen that H„ method 
provides a direct, systematic and effective way to synthesize controllers for the 
statically unstable longitudinal dynamics of X-29 aircraft. The two-command control 
design example demonstrates a fast response, large stability margins, good disturbance 
attenuation and low sensitivity performance. While the 3-input control configuration 
shows a somewhat less performance due to the smaller loop gains. 

It is seen in both design cases that the two control outputs, angle-of-attack a 
and pitch rate q, are decoupled to the command inputs. This is an important feature 
which is needed to effect the advanced control modes in which aircraft flight path and 
pitch attitude are independently controlled. Since, as shown in the 3-input case, the 
separated commanded strake has almost no influence on the control outputs, we may 
conclude that a reasonable and effective controller can be based on the use of two- 
command control in the X-29 control system . 

In both designs, we found the peak values of control surface deflections for the 
actuators and control rates are all much smaller than their design limits. The selected 
Wj weighting matrix ensures this compatibility but at a severe penalty in performance. 

Further reduced order modeling of the controller is possible but was not 
accomplished in this thesis. 
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APPENDIX A 



Matlab Script Files 



diary hx29.dat 
format short e 



disp( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp(' ') 
disp ( ' 
disp ( ' 
disp ( ' ') 

disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' ') 

disp ( ' ') 

%pause 
clc 
disp ( ' 
disp ( ' ') 

disp ( ' ' ) 

ag=[-. 4181d+00 .9960d+00 -.2269d-01 -.1213d+00 -.1948d-01 

. 4427d-04 -.6712d-04 .1451d-05 -.2162d-04 -.3540d-05 0.0 



') 

This script file is designed to solve the Hinf optimal control ') 
problem for the X-29. The 14th order FDLTI model, in state ') 
space form, is that of the X-29 aircraft and actuator dynamics.') 
Two states are those of the aircraft dynamics, i.e., alpha ') 
and q. The remaining 12 states are the dynamics of the ') 
three, fourth order actuators, i.e., the canard, flaperon, and ') 
strake actuators. The order of the unbalanced states is as ') 
follows : ' ) 

alpha, q, dc, df, ds, dcdot, dfdot, dsdot, dcdbldot,') 
dfdbldot, dsdbldot, dctrpldot, dftrpldot, dstrpldot') 

Given the open loop transfer function G ( s ) =Cinv ( Is-A) B+D , a ') 
stabilizing controller F(s) will be found such that the Hinf norm 
of Tylul is less than or equal to one. ') 



X-29 aircraft and actuator state space representation') 



• 54 74d + 0 1 -. 34 2 4d+00 



. 2585d+01 -. 1386d+01 



*. 1058d+01 



-.949 3 d-03 
0.0 0 . 0 ; 

. 3 898d-02 . 



. 1164d- 


-01 


-. 6397d-02 


- . 2509d-03 - 


. 53 62d-03 


-.2912d-< 


03 0. 


o 

o 

o 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 1.0 


0.0 


0 . 0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 . 0e+04 


0.0 0.0 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0e+04 0.0 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 . 0e+04 


0.0 

0.0 - 


0.0 - . 1479d+04 

. 2697d+03 0.0 


o 

o 

o 

o o 
o 


— . 


1143d+03 


o 

o 


0.0 


2529d+01 


0.0 . . 


0.0 

0.0 0 


0.0 
. 0 - 


0.0 - . 149 ld+04 0.0 

. 270 ld+03 0.0; 


0. 


0 -. 


. 1149d+03 


o 

o 


o 

o 


-.253 6d+0 1 .. 


0.0 0.0 
. 1790d + 02 


0.0 0.0 
0.0 0.0 - 


-. 5302d+05 
. 6053d+03 ] 


0. 


o 

o 


, 0 


18 16d+04 


0.0 


0.0 . 


• • 



bg=[0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



o; 

o; 

o? 

0 ; 

o; 

o; 



0 . 0 ; 
o.o; 

0 . 0 ; 

0 . 0 ; 

0 . 0 ; 

. 1479d+04 0.0; 

0.0 . 1491d+04 ; 

. 53 08d+05 0.0] 
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cg= [1.0 0.0 0.0 

0.0 1.0 0.0 



0.0 

0.0 



0.0 

0.0 



0.0 

0.0 



0.0 

0.0 



0.0 

0.0 



0.0 

0.0 



0.0 

0.0 



0.0 

0. 0 



0.0 

0.0 



0.0 

0.0 



0 . 0 ; 

0 . 0 ] 



dg=zeros (2) 

Ipause 

disp(' ') 

disp(' Balanced realization of the X-29 state space representation') 

disp(' ') 

[agbl , bgbl , cgbl, g , t ]=obalreal (ag , bg , eg) 
ag=agbl; bg=bgbl; cg=cgbl; 

Ipause 

w=logspace (-3,3, 100) ; 

svg=sigma (ag , bg , eg , dg , 1 , w) ; svg=20*logl0 (svg) ; 
semilogx(w, svg) 

title ('X-29 SV PLOT OF PLANT G') 
xlabel ( 'FREQUENCY - rad/sec') 
ylabel ( ' SV - dB') 
grid 

meta hx29gsv 
disp(' ') 

disp(' Calculate the poles and transmission zeros of the balanced open') 

disp(' loop plant') 

disp ( ' ') 

poleg=eig (ag) , tzerog=tzero (ag , bg , eg, dg) 
disp ( ' ') 

disp(' ') 

disp(' Determine determine the condition number of ag') 

disp(' ') 



disp ( ' ') 

condag=cond (ag) 



disp( ' 
disp ( ' 
Ipause 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp ( ' 



') 

') 



') 



1) 



rcondag=rcond (ag) 



<< Design Specifications >> ') 

Robustness Spec, : -40 dB roll-off, -20 db @ 100 Rad/Sec.') 
Associated Weighting:') 



-1 

W3(s) = 



2) . Performance Spec. 



1000 ') 

* I (f ixd) ' ) 

2 2x2') 

s ') 



minimizing the sensitivity function') 
as much as possible.') 

Associated Weighting:') 



') 



') 

-1 -1 . 01 (100s + 1) 

Wl(s) = Gam * * 

(.is + 1) 



') 

I') 

2X2') 

') 



where "Gam" in this design is iteratively updated from 1') 
w=logspace (-3 ,3,100) ; 
k=1000; mn= [ 2 2]; tau=0.0; 
nuv3i = [0.0 k] ; dnw3i = [1.0 0 0]; 
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sw3i = bode (nuw3i,dnw3i,w) ; svw3i = 2 0*logl0 ( svw3 i ) ; 
nuwli = [1.0 0.01]; dnwii =[0.1 1.0]; 

swli = bode (nuwli , dnwii ,w) ; svwii = 20*logl0 (svwli) ; 
aw2=-0.025*eye(2) ; bw2=zeros(2) ; cw2=zeros(2) ; dw2=-0 . 02 5* eye ( 2 ) ; 
disp(' ') 
disp ( ' ' ) 

disp(' (strike a key to see the plot of the weightings ...)') 

%pause 

semilogx (w, svwli , w, svw3i) 
grid 

title ('X-29 Design Specifications') 
xlabel (' Frequency - rad/sec') 
ylabel ( ' 1/W1 & 1/W3 - db') 
text ( .01,0, ' 1 /Wl ( s ) ') 



and find a stabilizing controller 
of TF Tylul is less than or equal 



) 



where ' ) 



text ( .5,100, ' 1/W3 (s) 
meta hx29specs 
%pause 
clc 
disp ( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp ( 9 
disp ( 9 
disp ( 9 
disp ( 9 
disp ( 9 
disp ( 9 
disp ( ' 
disp ( 9 
disp ( 9 
disp( ' 
disp ( 9 
Ipause 
clc 
disp ( 9 
disp ( 9 
disp( ' 
disp ( 9 
disp ( 9 
disp ( 9 
disp ( 9 
disp ( 9 
disp ( 9 
disp( ' 
disp ( 9 
disp ( 9 
disp ( ' 
disp( ' 
disp ( 9 
disp ( 9 
disp ( 9 



) 



<< Problem Formulation >>') 

Form an augmented plant P(s) with these two weighting functions:') 

1) . Gam*Wi penalizing error signal "e”') 

2) . W3 penalizing plant output M y M ') 

the Hinf norm') 



F ( s ) such that 
to one , i . e . ' ) 



F ( s) 



| Tylul | 



< or 
inf' ) 



= 1 ,') 



Tylul = 



') 

') 

') 



') 



') 



Gam*Wi* (I 
W3*GF* (I 



-1 

GF) 

-1 

GF) 



') 



') 



Gam 

W3 



Wl * 

(I - 



S 

S) 



') 

') 



(strike a key to continue ...)') 



« DESIGN PROCEDURE »') 

************************** 

[Step 1]. Do plant augmentation (run AUGMENT. M or 
AUGX29.M) 

(Step 2]. Balance the augmented plant for better 
numerical condition if necessary 

3] . Do Hinf synthesis with "Gam" = 1 

4] . Redo the plant augmentation for a 
higher M Gam" and rerun HINF.M 



[Step 

(Step 



********************* 



****** 



*') 

*') 

*') 

*') 

*') 

*') 

*') 

*') 

*') 

*') 

*') 

*') 
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disp ( ' 
disp ( ' 
%pause 
clc 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
haugx29 
disp ( ' 
disp ( ' 
disp ( ' 
%pause 
clc 
disp ( ' 
disp ( # 
%disp ( ' 
%disp ( ' 
%disp ( ' 



') 



') 

') 

') 



(strike a key to continue ...)') 



Assign the cost coefficients "Gam" with Gam=l 
serving as the baseline design 



augment 



% Plant augmentation of the X-29 dynamics') 



(strike a key to continue 



) ') 



) 

') 



Do state space balancing on the augmented plant if needed') 



), C2=cbal(5:6 / : ) 



% [abal , bbal , cbal , g , t ] =obalreal (A, [B1 B2] , [C1;C2]) 

%A=aba 1 , Bl«bbal(: , 1:2) , B2=bbal ( : , 3 : 4 ) , Cl=cbal(l:4 / 

') 

') 

The transmission zeros, poles and condition number of 
plant follow. In addition, determine if (A,B1) & (A,B2) 
stabilizable and if (Cl, A) & (C2,A) are detectable.') 

') 



’) 



disp( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 

tzeroaug=t zero ( A , [B1 B2 ] , [Cl ;C2 ] , [Dll D12;D21 D22]), poleaugA=eig (A) 
condaugA=cond (A) , rcondaugA=rcond (A) 
eps=eps 

toldef =10*max (size (A) ) *norm(A, 1) *eps 
tol=100*eps*norm ( [A Bl]) 

[Alc,Blc, Clc, t ,k]=ctrbf (A, Bl , Cl , tol ) 
tol=100*eps*norm ( [A B2]) 

[A2c,B2c,C2c, t,k]=ctrbf (A, B2 , C2 , tol) 
tol=100*eps*norm ( [ A; Cl ] ) 

[Alo,Blo,Clo,t,k] =obsvf (A, Bl , Cl , tol) 
tol=100*eps*norm ( [A;C2 ] ) 

[A2o, B2o , C2o, t , k]=obsvf (A, B2 , C2 , tol) 
clear condag rcondag poleg tzerog svw3i svwli 
clear condaugA rcondaugA 
clear functions 
%pause 

disp ( ' ') 

disp(' ') 

disp ( ' 

disp(' hinf 

disp ( ' 

aretype =/ Schur ' 
hinf 

disp(' ') 



the augmented') 
are ' ) 



% Running script file HINF . M for Hinf optimization') 
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(strike a key to continue ...)') 



disp(' ') 
disp( ' 

clear functions 
%pause 



disp ( ' 
disp( ' 
disp( ' 
disp ( ' 
disp( ' 



') 

') 

State 
(acp, 
) 



space representation of the full order controller') 
bcp, ccp # dcp) with its poles and condition number') 



acp , bcp , ccp , dcp 
polec=eig (acp) 
condacp=cond (acp) 
disp(' ') 

disp(' Minimal realization of the controller') 
disp(' ') 

toldef =10*max (size (acp) ) * norm (acp, 1) *eps 
tol=100*eps*norm ( [ acp bcp;ccp dcp]) 

[ acpm, bcpm , ccpm , dcpm] -minreal (acp , bcp , ccp , dcp) 
disp(' ') 

disp(' Balanced realization & model reduction of the minimal controller') 
disp(' ') 

[ acpbl , bcpbl , ccpbl , g , t ] =obalreal (acpm, bcpm , ccpm) 
elim= [ 10 ] 

[ acpr , bcpr , ccpr , dcpr ] =modred (acpbl , bcpbl , ccpbl , dcpm , elim) 
disp ( ' ') 

disp(' Poles, controllability, observability, and condition of the ') 

disp(' balanced, reduced order controller') 

disp ( ' ') 

poleacpr=eig (acpr) 

tol=100*eps*norm ( [acpr bcpr]) 

[acpc, bcpc , ccpc, t , k]=ctrbf (acpr, bcpr, ccpr , tol) 
tol=100*eps*norm ( [acpr ; ccpr ] ) 

[acpo, bcpo, ccpo, t , k ] =obsvf (acpr , bcpr , ccpr , tol ) 
condacpr=cond (acpr ) , rcondacpr=rcond (acpr) 
acp=acpr; bcp=bcpr; ccp=ccpr; dcp=dcpr; 
disp(' ') 

disp(' CLTF Tylul (acl, bcl, ccl, del) and its poles (reduced order)') 

[ acl , bcl , ccl , del ]=lf tf (sysp, dimp, acp, bcp, ccp, dcp) 
polet=eig (acl ) 

%pause 
disp ( ' ') 

disp ( ' ') 

hpltopmod 
end 

disp ( ' ' ) 

disp ( ' ' ) 

disp(' (strike a key to continue ...)') 

Ipause 
disp(' ') 
disp(' ') 
disp ( ' 
disp ( ' 



% Preparing singular values for plotting 



Open loop state space representation of controller/plant series') 
r ) 

[algf , blgf , clgf , dlgf ] = 
polol=eig (algf ) 



= series (acp, bcp, ccp, dcp, ag, bg, eg , dg) 



disp( ' 
disp( ' 
disp( ' 
%pause 
disp( ' 
disp( ' 



') 

') 



') 



(strike a key to continue ...)') 
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diary h3x29.dat 
format short e 



disp ( ' 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp ( 7 
disp( 7 
disp( 7 
disp ( 7 
Ipause 
clc 
disp( 7 
disp ( 7 
disp ( 7 



h3x29 . m 7 ) 

') 

This script file is designed to solve the Hinf optimal control 7 ) 
problem for the X-29. The 14th order FDLTI model, in state ') 
space form, is that of the X-29 aircraft and actuator dynamics. 7 ) 
Two states are those of the aircraft dynamics, i.e., alpha 7 ) 
and q. The remaining 12 states are the dynamics of the 7 ) 
three, fourth order actuators, i.e., the canard, flaperon, and 7 ) 
strake actuators. The order of the unbalanced states is as 7 ) 
follows : 7 ) 

alpha, q, dc, df, ds, dcdot, dfdot, dsdot, dcdbldot, 7 ) 
dfdbldot, dsdbldot, dctrpldot, dftrpldot, dstrpldot 7 ) 

Given the open loop transfer function G ( s) =Cinv ( Is-A) B+D, a 7 ) 
stabilizing controller F(s) will be found such that the Hinf norm 7 ) 
of Tylul is less than or equal to one. 7 ) 



') 



) 



X-29 aircraft and actuator state space representation 7 ) 



') 



r ) 



ag= [ - . 4 18 ld+00 . 9960d + 00 - . 22 69d-0 1 -.1213d+00 -.1948d-01 

.4427d-04 - . 67 12d-04 .1451d-05 -.2162d-04 -.3540d-05 0.0 



.5474 d+0 1 -.34 2 4d + 00 



. 2585d+01 - . 138 6d + 01 



*. 1058d + 01 



-. 9493d-03 
0.0 0 . 0 ; 

. 3898d-02 . 



1164d-01 


- . 6397d-02 


- . 2509d- 


•03 - 


. 5362d-03 


- . 2912d-i 


03 0. 


0 0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 . 0e+04 


0.0 0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0e+04 0.0; 


0.0 


0.0 


0.0 0.0 


0.0 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 . 0e+04 ; 


0.0 
. 0 - 


0.0 -.147 9d+04 

. 2697d+03 0.0 


0.0 0. 
0.0; 


0 -. 


1 14 3d+03 


o 

o 


o 

o 


- . 2 529d+01 


0.0 . . . 


0.0 
.0 0 


0.0 
. 0 - 


0.0 - . 1491d+04 0. 

■ . 2 701d+03 0.0; 


o 

o 


0 


114 9d + 03 


o 

o 


o 

o 


- . 2 536d + 01 . . . 


0.0 0.0 
1790d+02 


0.0 0.0 
0.0 0.0 - 


- . 5302d+05 0. 

. 6053d+03 ] 


0 0. 


0 -. 


18 16d + 04 


o 

o 


0.0 . 





bg= [0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 



0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 



0.0 0.0 
0.0 0.0 
0.0 0.0 
0.0 0.0 
0.0 0.0 
0.0 0.0 
. 14 79d+04 0.0 

0.0 .149 ld+04 



0 . 0 ; 

o.o; 



0.0 0.0 . 53 08d+05 ] ; 
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cg= [1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0 . 0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0 . 0 


0.0 


0.0 


0.0]; 



dg=zeros (3,3) ; 



%pause 

disp(' ') 

disp(' Balanced realization of the X-29 state space representation') 

disp(' ') 

[agbl , bgbl , cgbl , g, t]=obalreal (ag, bg, eg) 
ag=agbl , bg=bgbl, cg=cgbl, 
w=logspace (-3 ,3,100) ; 

svg=sigma (ag , bg , eg , dg , l,w) ; svg=20*logl0 (svg) ; 
semilogx (w, svg) 

title ('X-29 SV PLOT OF PLANT G') 
xlabel ( 'FREQUENCY - rad/sec') 
y label ( 'SV - dB' ) 
grid 

%meta h3x29gsv 
disp ( ' ') 

disp(' Calculate the poles and transmission zeros of the balanced open') 

disp(' loop plant') 

disp ( ' ' ) 

poleg=eig (ag) , tzerog=tzero (ag, bg, eg, dg) 
disp(' ') 

disp(' ') 

disp(' Determine determine the condition number of ag') 

disp ( ' ' ) 

disp(' ') 

condag=cond (ag) 



disp ( ' 
disp ( ' 
%pause 
disp(' ') 
disp ( ' 
disp(' ') 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp(' ') 
disp ( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp ( ' 



1 ) 



') 



) 



2 ) 



rcondag=rcond ( ag) 



<< Design Specifications >> ') 

Robustness Spec. : -40 dB roll-off, -20 db 0 100 Rad/Sec.') 
Associated Weighting:') 

-1 1000 ') 

W3 (s) = * I (f ixd) ' ) 

2 3x3') 

s ') 



Performance Spec.: minimizing the sensitivity function') 
as much as possible.') 

Associated Weighting:') 



') 

-1 -1 . 01 ( 100s + 1) 

Wl(s) = Gam * * 

( . Is + 1) 



I') 

3x3 ' ) 

') 



where ,, Gam H in this design is iteratively updated from 1') 



w=logspace (-3 , 3 , 100) ; 
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k=lOOO; mn=[3 3]; tau=0.0; 

nuw3 i = [0 0 k] ; dnw3i = [1.0 0 0]; 

svw3i = bode (nuw3i , dnv3i , w) ; svw3i = 20* loglO (svw3i) ; 
nuwii = [l.o 0.01]; dnwli =[0.1 1.0]; 

svwli = bode (nuwii , dnwli ,w) ; swli = 20*logl0(svwli); 
aw2=-0. 025*eye (3) ; bw2=zeros(3) ; cw2=zeros(3) ; dw2=-0 . 02 5*eye ( 3 ) ; 
disp ( ' ') 

disp(' ' ) 

disp(' (strike a key to see the plot of the weightings ...)'; 

%pause 

semilogx(w, svwli, w # svw3i) 
grid 

title ( ' X-2 9 Design Specifications') 
xlabel (' Frequency - rad/sec') 
ylabel ( ' l/Wi & 1/W3 - db') 
text (.01,0, '1/Wl(s) ') 
text ( . 5 , 100 , ' 1 /W3 ( s ) ') 

%meta h3x29spc 

%pause 

clc 

disp(' << Problem Formulation >>') 

disp ( ' ') 

disp(' Form an augmented plant P(s) with these two weighting functions:') 
disp ( ' ' ) 

1) . Gam*Wl penalizing error signal "e"') 

') 

2) . W2 penalizing control input "u"') 

') 

2). W3 penalizing plant output "y"') 

') 

and find a stabilizing controller F(s) such that the Hinf norm') 
of TF Tylul is less than or equal to one, i.e.') 



disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' ' 

disp ( ' 
disp ( ' 
disp(' ' 
disp(' ' 
disp ( ' 
%pause 
clc 

disp(' ' 
disp(' ') 
disp ( ' 
disp(' ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp( ' 
disp(' ' 
disp ( ' ' 

disp ( ' 
Ipause 
clc 



| Tylul | < or 

F ( s) inf') 



1 ,') 



(strike a key to continue ...)') 



<< DESIGN PROCEDURE >>') 

******************************') 

* [Step 1]. Do plant augmentation (run H3AUGX29.M ) * ' ) 



★ 

* 

* 

* 

★ 

★ 

* 

* 

******* 



[Step 2]. Balance the augmented plant for better 
numerical condition if necessary 

[Step 3]. Do Hinf synthesis with "Gam” = 1 

[Step 4]. Redo the plant augmentation for a 
higher n Gam M and rerun HINF.M 



*') 

* # ) 

*') 

*') 

*') 

*') 

*') 



**********************') 

(strike a key to continue ...)') 
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disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( 9 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 
h3augx29 
disp ( ' 
disp ( ' 
disp ( 9 
%pause 
clc 
disp ( 9 
disp( ' 
%disp ( ' 
%disp ( 9 
%disp ( 9 



Assign the cost coefficients "Gam" with Gam=l 
serving as the baseline design 



augment 



% Plant augmentation of the X-29 dynamics') 



(strike a key to continue 



) ') 



Do state space balancing on the augmented plant if needed') 



: ) , C2=cbal(5:6, : ) 



% [abal , bbal , cbal , g, t]=obalreal (A f [B1 B2 ] , [C1;C2] ) 

%A=abal , Bl=bbal ( : , 1 : 2 ) , B2=bbal ( : , 3 : 4 ) , Cl=cbal(l:4, 

') 

') 

The transmission zeros, poles and condition number of the augmented') 
plant follow. In addition, determine if (A,B1) & (A,B2) are ') 
stabilizable and if (Cl, A) & (C2,A) are detectable.') 



disp ( ' 
disp( ' 
disp ( ' 
disp ( ' 
disp ( ' 
disp(' ') 
disp(' ') 

tzeroaug=tzero (A, [B1 B2 ] , [Cl ; C2 ] , [Dll D12;D21 D22]), poleaugA=eig (A) 

condaugA=cond (A) , rcondaugA=rcond (A) 

eps=eps 

toldef =10*max (size (A) ) * norm (A, 1) *eps 
tol=100*eps*norm ( [A Bl]) 

[Alc,Blc,Clc,t,k] =ctrbf (A , Bl , Cl , tol ) 
tol=100*eps*norm( [A B2]) 

[ A2c , B2c , C2c , t , k] =ctrbf ( A, B2 , C2 , tol ) 
tol=100*eps*norm( [A; Cl] ) 

[Alo,Blo,C10,t,k] =obsvf (A, Bl, Cl, tol) 
tol=100*eps*norm ( [A; C2 ] ) 

[A2o, B2o, C2o, t,k]=obsvf (A, B2 , C2 , tol) 
clear condag rcondag poleg tzerog svw3i svwli 
clear condaugA rcondaugA 
clear functions 



%pause 



disp(' ') 

disp(' ') 

disp ( ' ') 

disp(' hinf % Running script file HINF.M for Hinf optimization') 

disp ( ' ') 



aretype=' Schur ' 
hinf 

disp( ' ' ) 

disp(' ') 

disp(' (strike a key to continue ...)') 

clear functions 

%pause 
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) 



space representation of 
bcp, cep, dep) with its 



the full order controller') 
poles and condition number') 



and condition of the ') 



disp ( ' 
disp(' ') 
disp ( ' State 

disp(' (acp, 

disp ( ' ' ) 

acp, bcp, ccp, dep 
polec=eig (acp) 
condacp=cond (acp) 
disp(' ') 

disp(' Minimal realization of the controller') 
disp(' ') 

toldef =10*max (size (acp) ) *norm (acp, 1) *eps 
tol-100*eps*norm ( [ acp bcp;ccp dep]) 

[acpm, bepm, ccpm, depm] =minreal (acp, bcp, ccp , dep) 
disp ( ' ') 

disp(' Balanced realization & model reduction of the minimal controller') 
disp ( ' ') 

[acpbl,bcpbl,ccpbl,g, t]=obalreal (acpm, bepm, ccpm) 
elim= [ 10 ] 

[aepr , bepr , ccpr , depr ]=modred (acpbl , bcpbl , ccpbl , depm, elim) 
disp(' ') 

disp(' Poles, controllability, observability, 

disp(' balanced, reduced order controller') 

disp(' ') 

poleacpr=eig (aepr) 
tol=100*eps*norm ( [ aepr bepr]) 

[aepe, bepe, ccpc, t,k]=ctrbf (aepr, bepr, ccpr, tol) 
tol=100*eps*norm ( [aepr ; ccpr ] ) 

[aepo, bepo, ccpo, t,k]=obsvf (aepr, bepr, ccpr, tol) 
condacpr=cond (aepr) , rcondacpr=rcond (aepr ) 
acp=acpr; bcp=bcpr; ccp=ccpr; dcp=dcpr; 
disp(' ') 

disp(' CLTF Tylul (acl, bcl, eel, del) and its poles (reduced order)') 
[acl,bcl,ccl,dcl]=lftf (sysp, dimp, acp, bcp, ccp, dep) 
polet=eig (acl) 

%pause 
disp( ' 
disp ( ' 
h3plot 
end 
disp ( ' 
disp ( ' 
disp ( ' 

%pause 
disp ( ' 
disp ( ' 
disp ( ' 
disp ( ' 



') 

') 



') 



') 

') 



% Preparing singular values for plotting 



(strike a key to continue ...)') 



Open loop state space representation of controller/plant series') 
) 

[ algf , blgf , clgf , dlgf ] = series (acp, bcp, ccp, dep, ag, bg, eg, dg) 
polol=eig (algf) 



disp ( ' 
disp ( ' 
disp ( ' 
Ipause 
disp ( ' 
disp ( ' 
disp ( ' 
disp( ' 
disp ( ' 
disp( ' 



') 

') 



’) 



(strike a key to continue 



) ') 



') 

Closed loop state space representation of controller/plant series,') 
controllability, observability, and condition number of the closed') 
loop aegf matrix. ') 

') 
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disp ( ' ' ) 

[acgf , bcgf , ccgf , dcgf ] =f eedbk (algf , blgf,clgf,dlgf,2) 
tol=100*eps*norm ( [acgf bcgf]) 

[acgf c, bcgf c, ccgf c, t , k] =ctrbf (acgf , bcgf , ccgf , tol ) 
tol=100*eps*norm( [acgf ; ccgf ] ) 

[acgf o, bcgf o, ccgf o, t , k] =obsvf (acgf , bcgf , ccgf , tol) 
condacgf=cond (acgf ) 
disp(' ') 
disp(' ') 

disp(' (strike a key to continue ...)') 

%pause 

disp ( ' ') 

disp(' Poles of the closed loop system') 

disp(' ') 

polcl=eig (acgf ) 

h3aly 

%h3rsp 

end 
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% 

% Plant Augmentation for the X-29 Hinf problem as W3 is not a 
% proper transfer function. Includes contingency for adding W2 to 
% ensure dl2 is full column rank. This script file is designed for 
% the X-29 system without theta as a state, ie, 3 inputs & 3 outputs. 

%disp ( ' ') 

%disp(' ') 

%disp(' << Plant Augmentation >>') 

%Gam=gama ( 1 , i) 

Gam = input (' Input the cost coefficient "Gam" = ') 

cgb = 1/k* [ eg ( 1 , : ) *ag*ag ; 

eg (2 , : ) *ag*ag*ag*tau+cg (2,:)*ag*ag; 
eg (3 , : ) *ag*ag*ag*tau+cg (3 , : ) *ag*ag ] ; 
dgb = 1/k* [ eg ( 1 , : ) *ag*bg ; 

eg (2 , : ) *ag*ag*bg*tau; 
eg (3 , : ) *ag*ag*bg*tau] ; 

nwl = Gam*[dnwii;0 0;0 0;0 0;dnvii;0 0;0 0;0 0;dnwli]; 
dwl = nuwli; 

sysw2=[aw2 bw2;cw2 dw2]; xw2=3; 

%[A,B1, B2 ,C1,C2 ,D11,D12 ,D21,D22] =haugmod (ag , bg , eg , cgb, dg , dgb , nwl , dwl , mn) 

[ A, B1 , B2 , Cl , C2 , Dll , D12 , D21 , D2 2 ] =haugx2 9 pi (ag, bg , eg , cgb , dg , dgb , nwl , dwl , sysw2 , xw2 , 

%disp ( ' ') 

%disp ( ' - - - State-Space (A, B1,B2 , Cl, C2 , Dll, D12 , D21, D22) is ready for') 

%disp(' the Small-Gain problem - - -') 

% 

% End of AUGX29.M % A Z 
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disp(' Closed loop state space representation of controller/plant series,') 
disp(' controllability, observability, and condition number of the closed') 
disp(' loop acgf matrix. ') 
disp(' ') 
disp(' ') 

[acgf , bcgf , ccgf , deg f]=feedbk (algf , blgf , clgf,dlgf,2) 
tol=100*eps*norm( [acgf bcgf]) 

[acgfc, bcgf c, ccgf c,t,k]=ctrbf (acgf , bcgf , ccgf, tol) 
tol=100*eps*norm ( [ acgf ; ccgf ] ) 

[acgf o, bcgf o, ccgf o, t , k]=obsvf (acgf, bcgf , ccgf , tol) 
condacgf =cond (acgf) 
disp(' ') 
disp ( ' ' ) 

disp(' (strike a key to continue ...)') 

Ipause 

%diary hx29.dat 
disp(' ') 

disp(' Poles of the closed loop system') 
disp(' ') 
polcl=eig (acgf ) 

%hx29aly 

%hx29trsp 

end 
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format short e 
w=logspace ( -3 , 3 , 100) ; 

[magi , phasel ]=bode (algf , blgf , clgf , dlgf , l,v) ; magl=20* log 10 (magi) 

[mag2 , phase2 ] =bode ( algf , blgf , clgf , dlgf , 2 , w) ; mag2=20*logl0 (mag2 ) 

[mag 3 , phase3 ] =bode (algf , blgf , clgf , dlgf , 3 , v) ; mag 3=2 0* log 10 (mag3 ) 

semilogx (v, magi ( : , 1) , w,mag2 ( : , 1) , w, mag3 ( : , 1) ) 

title ( ' X-29 OPEN LOOP (GF) BODE PLOT OF " alpha " ') 

xlabel(' FREQUENCY - rad/sec') 

ylabel ( ' GAIN - dB') 

grid 

%meta h3bdl 
%pause 

semi logx(w, magi ( : , 2) ,w,mag2 ( : , 2) ,w / mag3 ( : , 2) ) 
title (' X-29 OPEN LOOP (GF) BODE PLOT OF M q " ') 
xlabel( 'FREQUENCY - rad/sec') 
ylabel ( 'GAIN - dB' ) 
grid 

%meta h3bd2 
%pause 

[magi , phasel j =bode (acgf , bcgf , ccgf , dcgf , 1 , w) ; magl=20*logl0 (magi) 

[mag2 , phase2 ] =bode (acgf , bcgf , ccgf , dcgf , 2 , w) ; mag2=20*logl0 (mag2 ) 

[mag3 , phase3 ]=bode ( acgf , bcgf , ccgf , dcgf , 3 , w) ; mag3=20*logl0 (mag 3 ) 

semilogx ( w, magi ( : , 1) , w ; mag2 ( : , 1) , w ; mag3 ( : , 1) ) 

title ('X-29 CLOSED LOOP BODE PLOT OF " alpha M ') 

xlabel( 'FREQUENCY - rad/sec') 

ylabel ( ' GAIN - dB') 

grid 

text (15, -20, ' TO rl') 

text ( .22,-20, ' TO r2') 

text (.1,-60, ' TO r3 ' ) 

meta h3bd3 
pause 

semilogx (w, magi ( : ,2) ,w,mag2 ( : ,2) ,v,mag3 ( : ,2) ) 
title ('X-29 CLOSED LOOP BODE PLOT OF " q " ' ) 
xlabel ( 'FREQUENCY - rad/sec') 
ylabel ( 'GAIN - dB' ) 
grid 



text(. 1 ,- 10 , 


--TO 


r2 ' ) 


text(. 1,-30, ' — 


TO 


rl') 


text ( . 1 , -58 , ' — 


TO 


r3') 



meta h3bd4 
%pause 
disp ( ' ') 

end 
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APPENDIX B 



State Space Realization of the Open Loop Plant G 



ag=[-. 4181d+00 .9960d+00 -.2269d-01 -.1213d+00 -.1948d-01 -.9493d-03 .. 



4427d- 


04 - 


. 67 12d-04 


. 14 51d-05 




2 1 62d- 


04 - 


. 3 54 0d-0 5 


0.0 


0.0 


0.0; 


. 5474d+01 




34 24d + 00 


. 2 58 5d+0 1 


-. 138 6d + 01 


1058d+01 


. 3898d-02 ... 


. 1164d 


-01 


6397d-02 




2 509d- 


03 - 


-. 5362d-03 


2912d-i 


03 0. 


0 0.0 


o.o; 


0.0 


0.0 


0. 0 


0.0 


0. 0 


1.0 


0.0 


0.0 


0.0 


0. 0 


0.0 


0 . 0 


0.0 


0.0; 


0.0 


0.0 


0.0 


0.0 


0. 0 


0.0 


1 . 0 


0.0 


0.0 


0.0 


0.0 


0. 0 


0.0 


0.0; 


0.0 


0.0 


0. 0 


0. 0 


0.0 


0.0 


0 . 0 


1 . 0 


0. 0 


0.0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0. 0 


0. 0 


0.0 


0.0 


0.0 


0 . 0 


0.0 


1 . 0 


0.0 


0.0 


0. 0 


0.0 


0.0; 


0. 0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 . 0 


0.0 


0.0 


0.0 


0.0; 


0.0 


0.0 


0.0 


0. 0 


0.0 


0.0 


0. 0 


0.0 


0.0 


0. 0 


1.0 


0.0 


0.0 


0.0; 


0.0 


0. 0 


0.0 


0.0 


0.0 


0.0 


0. 0 


0.0 


0.0 


0. 0 


0.0 


1 . oe+04 


0.0 0.0 


0.0 


0.0 


0.0 


0. 0 


0.0 


0.0 


0. 0 


0.0 


0.0 


0. 0 


0.0 


0 . 0 


1.0e+04 0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 . 0e+04 


0.0 


0.0 


14 7 9d + 04 


0 


. 0 0. 


0 - 


. 1143d+03 


0.0 


0. 0 


2529d+01 


0.0 . . 


o.o - 


. 2697d + 03 


0.0 


0.0 


/ 


















0.0 


0.0 


0.0 


- . 1491d+04 0. 


0 0 


. 0 -. 


114 9d 


+ 03 


0.0 


0.0 


- . 2 53 6d+01 .. 


0.0 0 


. 0 - 


. 2701d + 03 


0.0 


t 


















0.0 


0.0 


0.0 


0.0 


- . 53 02d + 05 0 


. 0 0. 


0 


181 6 d +04 


0. 0 


0.0 . 


• • 



- . l790d + 02 0.0 0.0 -,6053d+03) 



bg= [0.0 0.0; 

0.0 0 . 0 ; 
o.o o.o; 

0.0 0 . 0 ; 

O.o o.o; 

0.0 0.0; 

0.0 0.0; 

0.0 0.0; 

0.0 0.0; 

0.0 0.0; 

0.0 0.0; 

. 1479d+04 0.0; 

0.0 . 149 ld+04 ; 

. 53 08d+05 0.0] 



cg= [1.0 0.0 

0.0 1.0 

dg=zeros (2 ) 



0.0 0.0 0.0 

0.0 0.0 0.0 



0.0 0.0 0.0 

0.0 0.0 0.0 



0.0 0.0 0.0 

0.0 0.0 0.0 



0.0 0.0 0.0; 

0.0 0.0 0 . 0 ] 
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APPENDIX C 



State Space Realization of Reduced-Order Controller for 2 I/O System 



acpr = 



Columns 1 through 6 



-1. 0399e-02 
-1. 127le-04 
-1.0340e-01 
-5.7615e-02 
7. 8601e-02 
-3 . 4 04 8e-03 



-1 . 22 62e-05 
-1 . 0156G-02 
-2 . 9977e-02 
6 . 0104e-02 
1. 3419e-02 
-3 . 61 38e-02 



-1. 0460e-01 
**4 . 679 5e-03 
-1.9525e+01 
-2.7029e+01 
3 . 1746e+01 
-2 . 7492e+00 



4 . 8094G-02 
8 . 3252G-02 
2 . 6394e+01 
-2 . 5674e+01 
3 . 64 83e+01 
5 . 9997e+00 



-7 .68 7 5e-02 
-3 . 18 61e-02 
-3. 1280e+01 
4 . 0396e+01 
-1 . 2030e+02 
-2 . 9572e+01 



1. 3274e-02 
-2 . 8000G-02 
4 . 92 71e+00 
-7.7033e-02 
5. 3509e+01 
-9.8984e+00 



-4 . 98 64e-03 


-3. 9735e-02 


-3. 1302e+00 


3 . 5764e+00 


5 . 9 2 53e+01 


-2.4 7 77e + 01 


3 . 6535e-03 


-2. 9112e-02 


2 . 9037e-01 


1 . 0513e+01 


-2 . 47 17e+01 


-2 . 1 68 5e + 01 


4.79 08e-02 


-1. 1300G-02 


1 . 7010e+01 


3.28 1 3e+01 


-1.4 766e+02 


1 . 029 0e + 01 


9 . 2 037e-03 


-1 . 4 652e-03 


3 . 2877e+00 


6. 3758e+00 


-2 . 8661e + 01 


8 ♦ 9049e-01 


-8 . 7 2 00e-04 


6 . 8844e-03 


-9 . 962 le-02 


-2. 1329e+00 


1 . 64 4 6e + 00 


3 . 5846e+00 


8 . 304 6e-04 


-3 . 7921e-03 


1 . 8 I96e-01 


1 . 4 002e+00 


-1.9663e+00 


-1 . 8969e+00 


8 . 244 5e-04 


1 . 094 2e-04 


3 . 0302e-01 


5. 1368e-01 


-2 . 52 18e+00 


1 . 7864e-01 


-1.7763e-04 


9 . 0089G-04 


-3. 6094G-02 


-3 . 198 le-01 


4 . 0668e-01 


4 . 54 02e-01 


2 . 3 188e-05 


3.6527G-05 


9. 5788G-03 


6.88 3 2e-03 


-7 . 5896e-02 


2. 2842e-02 


Columns 7 through 12 










2.5870e-03 


9.1286G-03 


8 . 38 8 6g- 03 


8 . 694 0e-03 


5. 3095e-04 


1 . 7503e-03 


-4.0336G-02 


2 . 5238e-02 


-8.4033G-02 


-5. 5085e-03 


-7 . 0034e-03 


-2 . 6531e-03 


9.2118e-0l 


3 . 4965e+00 


2 . 3161e+00 


3. 1564e+00 


1. 5463e-01 


6.28 54 e-01 


1 . 014 2e+01 


-1 . 3763e+01 


1.9246e+01 


-3. 0766e+00 


1 . 85l0e+00 


-1. 3642e-01 


-5. 1779e+01 


4 .6467e+01 


1.98 5 3e+0 1 


2 . 4744e+01 


-1 . 1656e+00 


4 . 1721e+00 


3.2298G+00 


1 . 4770e+01 


-4 . 2332e+01 


-7. 6418e+00 


-3 . 28 32e+00 


-2 ♦ 1705e+00 


-1.3998e+01 


3.5254e+01 


-6 . 3778e+01 


-5 . 0628e+00 


-4 . 9278e+00 


-2. 0558e+00 


-3.7012e+01 


-1.4396e+01 


9 . 7 679e + 01 


-3 . 02 19e+00 


5. 0051e+00 


5 . 6593e-01 


2 . 04 92e+01 


-7 . 64 98e+01 


-2 . 4 305e+02 


-2.9439e+0l 


-4 . 1233e+01 


-1.7944e+01 


1. 5648e+00 


-4 .6369e+00 


-8 . 7073e+01 


-4 . 5600e+01 


-2 . 4 619e+01 


-2. 1083e+01 


4 . 4792e+00 


6.6710e+00 


1. 6807e+01 


2 . 0142e+01 


-2 . 72 37e+01 


-2.58 9 5e+01 


-2 . 3356e+00 


-3 . 89 24e + 00 


-1. 1275e+01 


-9 . 79 36e + 00 


3 . 29 61e + 0 1 


-1 . 2550e+02 


2 . 7969e-01 


-2. 1520G-01 


-8. I078e+00 


-8. 0534e+00 


1 . 9334e+00 


-2 . 7568e+01 


5.6035G-01 


9.1670G-01 


2 ♦ 5024e+00 


2. 2396e+00 


-7. 3070e+00 


5 . 738 0e+01 


3.0310G-02 


2 . 54 82e-02 


-2 . 0105e-01 


-2. 0731e-01 


-2 . 0302e-01 


1. 1700e+00 


Columns 13 through 15 










-4.8787e-04 


-4 . 7600e-04 


-1. 1648e-05 








-1. 1798e-03 


-4 . 622 3e-04 


-5 . 0636e-05 








-1.864 le-01 


- 1 . 7 78 5g-0 1 


-4 . 58 15e-03 








6 . 2 1 64 g-01 


3.9705G-01 


2. 1676e-02 








-1. 9086e+00 


-1. 5949e+00 


-5.4 32 9e-02 








-1 . 92 15e-01 


9 . 88 20e-02 


-1. 3922e-02 








-7 ♦ 37 4 2e-01 


-2 . 4 88 6e-01 


-3.2950G-02 








1 . 4035e+00 


8.0952G-01 


5. 1750e-02 








-5 . 7313e+00 


-1.74 34e+00 


-2 . 6308e-01 








2 . 34 04e+00 


3 . 5305e+00 


1. 1476G-02 








-8.154 3e+00 


-2 . 5954e+00 


-3 . 651 le-01 








2.4647e+00 


1 . 5106e+01 


-1 . 0758e+00 








-4 .5631e+0l 


-6 . 34 4 0e + 01 


-3 ♦ 3 130e+00 








3. 3396e+01 


-7. 1315e+01 


-6. 2957e+00 








-2 . 7014e+00 


-3.0793e+00 


-1 . 44 18e+02 
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bcpr = 



5 . 2 502e + 00 


6. 4618e-01 


-3 . 2830e-01 


2 . 9101e+00 


2 . 5970e+01 


7.309 le+00 


1 . 5589e+01 


“7. 1255e+00 


-1.9 609e + 01 


-4 . 1407e+00 


2 . 1893e-01 


5. 2816e+00 


5 . 578 3e-01 


5 . 8 4 52 e + 00 


“1 . 4366e+00 


4 . 0563e+00 


-1.2 3 10e + 01 


2 . 5979e“01 


-2 . 3 52 2e+00 


-5.2653e-02 


3 . 4 165e-01 


-9 . 602 7e-01 


-2 . 7670e-01 


5. 1899e-01 


-2 . 064 6e“01 


-3 . 9 167e-02 


6. 0765e-02 


“1 . 2386e-01 


-5 . 2 177e-03 


-5.8883e-03 



ccpr = 




Columns 1 through 6 


“3 . 8779e+00 


1 . 982 9e+00 


“3 . 5979e+00 


-2 . 1551e+00 


Columns 7 through 12 


4 ♦ 4 501e+00 


-7 . 6160e-01 


“3. 8307e+00 


4 . 2 3 52e + 00 


Columns 13 ‘ 


through 15 


2 . 4 380e-02 


“4 . 3767e-02 


*2 . 0872e-0 1 


-I. 3083e-01 


dcpr = 




-3 . 9822e-01 


-9. 2951e-02 


-6. 3988e-02 


“1. 4936e-02 



-1 . 9257e+01 
-1. 8895e+01 



9 . 8134e+00 
-7 . 4360e+00 



2 . 7852e-03 
-7.3580e-03 



7 . 12 36e-01 
1. 7126e+01 



2 . 1675e+00 
9. 1514 e-01 



-1. 1251e+01 
-1. 6585e+01 



7 . 8 612e-01 
“6.4874e-01 



5.24 13e+00 
6 • 8 677e-01 



5 . 8776e-01 
2 . 1019e-02 
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APPENDIX D 



Closed Loop State Space Realization of Controller/Plant Series 

Columns 1 through 6 



-1. 0399e-02 


-1. 2262e-05 


-1 . 04 60e-01 


4 . 8094 e-02 


-7 . 687 5e-02 


1. 3274e-02 


-1. 127le-04 


-1. 0156e-02 


-4 . 6795e-03 


8 . 3 2 52e-02 


-3 . 1861e-02 


-2 . 8000e-02 


-1.03 4 0e-01 


-2 . 9977G-02 


-1. 9525e+0l 


2 . 63 9 4 e + 01 


-3 . 12 80e + 01 


4 .927 le+00 


-5 . 7 615e-02 


6 . 0104e-02 


-2.7 02 9e + 01 


-2 . 5674e+01 


4 . 0396e+01 


-7 . 7 03 3e-02 


7 . 8601e-02 


1. 34 19e-02 


3 . 174 6e + 01 


3 . 64 83e+01 


-1 . 2030e+02 


5 . 3 509e+01 


-3 . 4 04 8e-03 


-3 . 613 8e-02 


-2 . 7492e+00 


5.9997e+00 


-2 . 9572e+01 


-9 . 8984e+00 


—4 . 9864e-03 


-3 . 97 35e-02 


-3 . 13 02 e+00 


3 . 57 64e + 00 


5 . 92 53e + 01 


-2 . 4777e+01 


3 . 653 5e-03 


-2 . 9112e-02 


2 . 9037e-01 


1 . 0513e+01 


—2 . 47 17e+01 


-2 . 1685e+01 


4 .7908e-02 


-1. 13 00e-02 


1 . 7 0 10e+01 


3 . 2813e + 01 


-1. 4766e+02 


1.0290e+01 


9.2037e-03 


-1. 4652e-03 


3 . 2877e+00 


6 . 3758e+00 


-2 . 8661e+01 


8 . 904 9e-01 


-8 . 7 200e-04 


6. 8844G-03 


-9 . 962 le-02 


-2 . 1329e+00 


1. 64 4 6e + 00 


3 . 5846e+00 


8. 3046e-04 


-3.792 le-03 


1 . 8196e-01 


1 . 4 002e + 00 


-1 . 9663e+00 


-1 . 89 69e + 00 


8.2445e-04 


1 . 094 2e-04 


3 . 0302e-01 


5 . 13 68e-01 


-2 . 52 18e+00 


1 . 7864e-01 


-1.77 63e-04 


9 . 0089e-04 


-3 . 6094 e-02 


-3 . 1981e-01 


4 . 0668e-01 


4 . 54 02e-01 


2 . 3188e-05 


3 . 652 7e-05 


9 . 5788e-03 


6 . 8832e-03 


-7.5896e-02 


2.284 2 e-02 


-1 . 2486e+00 


-2 . 9277e-01 


-6.4 4 00e+00 


4 . 067 5e + 00 


-5 . 0553e+00 


3 . 9391e-01 


6. 2283e-01 


-3. 5136e+00 


2 . 27 04e + 00 


1. 3055e+01 


-3. 1083e+00 


-5 . 2805e+00 


-2 . 82 3 0e+00 


4 . 3606e+00 


-1. 32 68e+01 


-1.1504e+01 


-3 . 7 03 2e + 00 


7 . 8 67 9e+00 


4 . 5916e-01 


-1 . 82l7e+00 


1 . 87 16e+00 


6 . 4561e+00 


-1 . 1091e+00 


-2 . 8251e+00 


-2 . 6795e+00 


7 . 004 9e-01 


-1.3478e+01 


3 . 252 3e + 00 


-8 . 8046e+00 


2 . 6913e+00 


3 . 34 98e + 00 


-6 . 0600e-01 


1. 6920e+01 


-5. 1776e+00 


1. 1422e+01 


-2 . 9899e+00 


-1. 1953e+00 


1 . 4 17 2e+00 


-5 . 7 2 85e + 00 


-3 . 102 3e + 00 


-2 . 2283e+00 


2 . 7 3 53e+00 


-1. 9331e-0l 


9 . 004 3e-01 


-7. 5360G-01 


-3 . 2683e+00 


6 . 7229e-01 


1. 3748e+00 


1. 1620e+00 


-1.4 09 3g-01 


5. 8870e+00 


-2 . 0816e+00 


4 . 0687e+00 


-9 . 4 087e-01 


-9 . 4 959e-0 1 


5 . 537 2e-01 


—4 . 698 le+00 


-1 . 064 le-01 


-2. 6504e+00 


1.3781e+00 


-4 . 6009e-0l 


1 . 092 2e-02 


-2 . 34 25e+00 


1. 0092e+00 


-1 . 6801e+00 


3. 102 0e-01 


1.4474 e-01 


7.964 6e-02 


7. 5834e-0l 


-6 . 599 le-01 


6 . 563 6e-01 


1 . 7830e-02 


1. 3788e-02 


-9.7 163G-02 


4 . 5274G-02 


3 . 6887e-01 


-9 . 862 3e-02 


-1 . 4 382e-01 


1 . 363 0e-02 


-6.9743e-03 


6 . 7 686e-02 


-2 . 4 859e-03 


3.954 le-02 


-1 . 8429e-02 



Columns 7 through 12 



2 . 587 0e-03 


9. 1286e-03 


8 . 3 88 6e-03 


8. 694 0e-03 


5. 3095e-04 


1 . 7503e-03 


-4 . 0336e-02 


2 . 523 8e-02 


-8 . 4 03 3e-02 


-5. 5085e-03 


-7 . 0034e-03 


-2 . 653 le-03 


9. 2118e-0l 


3 . 4 965e + 00 


2 . 3 161e+00 


3 . 1564e+00 


1 . 54 63e-01 


6 . 2854e-01 


1 . 014 2e+0 1 


-1. 3763e+01 


1. 9246e+01 


-3 . 0766e+00 


1 . 8510e+00 


-1. 3642e-01 


-5. I779e+01 


4 . 6467e+01 


1. 9853e+01 


2 . 4744e+01 


-1. 1656e+00 


4 . 1721e+00 


3 . 2298e+00 


1. 4770e+01 


-4 . 2 332e+01 


-7.64 18e+00 


—3 . 2 832e + 00 


-2. 1705e+00 


-1 . 3998e+01 


3 . 52 54e+01 


-6. 3778e+01 


-5 . 062 8e+00 


-4 . 92 7 8e + 00 


-2 . 0558e+00 


-3 . 7012e+01 


-1.439 6e + 01 


9 . 7 679e+01 


-3 . 02 19e + 00 


5 . 0051e+00 


5. 6593e-01 


2 . 04 92e+01 


-7 . 6498e+01 


-2 . 4 305e+02 


-2 . 94 39e + 01 


-4 . 12 33e + 01 


-1 . 7944e+01 


1 . 5648e+00 


-4 . 6369e+00 


-8 . 707 3e+01 


-4 . 5600e + 01 


-2 . 4619e + 01 


-2 . 1083e+01 


4 . 4792e+00 


6 . 67 10e+00 


1. 6807e+01 


2 . 014 2e+0l 


-2 . 7237e + 01 


-2 . 5895e+01 


-2. 3356e+00 


-3 . 8924e+00 


-1. 1275e+01 


-9 . 793 6e+00 


3 . 2961e + 01 


-1 . 2 550e + 02 


2 . 79 69e-0 1 


-2. 1520e-01 


-8. 1078e+00 


-8 . 0534e+00 


1 . 93 34e+00 


-2 . 7568e+01 


5 . 603 5e-01 


9 . 167 Oe-0 1 


2 . 5024e+00 


2 . 2 396e + 00 


-7.3 07 0e+00 


5 . 7380e+01 


3 . 03 10e-02 


2 . 54 82e-02 


-2 . 0105e-01 


-2 . 07 3 le-0 1 


-2 . 0302e-01 


1. l700e+00 


-4 . 2260e-01 


9. 0677e-01 


-6 . 9611e-01 


4 . 4243e-01 


-6 . 8 13 6e-02 


6 . 7 025e-02 


-7 . 0809e+00 


4 . 0749e+00 


- 1 . 4806e+01 


-1. 2247e+00 


-1 . 2285e+00 


-5.137 5e-01 


9 . 05l6e+00 


-4 . 163 0e+00 


1.9222e+01 


2 . 3782e + 00 


1. 5786e+00 


8 . 107 3e-01 


-3. 6887e+00 


2 . 053 3e+00 


-7.73 2 6e + 00 


-6. 9198e-01 


-6 . 4 051e-01 


-2 . 7787e-01 


1 . 7 4 06e + 00 


3.0217e-01 


4 . 0080e+00 


1. 314 0e+00 


3 . 12 17e-01 


3 . 1824e-01 


-1. 6387e+00 


-7 . 114 5e-01 


-3 . 89 38e + 00 


-1. 5687e+00 


-2.972 le-01 


-3 . 624 5e-01 


2 . 977 6e+00 


-1. 2318e+00 


6 . 3 62 2e+00 


8 . 8925e-01 


5 . 2036e-01 


2 . 8696e-01 


1 . 8 190e+00 


-1 . 029 6e + 00 


3 . 8082e+00 


3 . 2795e-01 


3 . 1571e-01 


1 . 3450e-01 


-4 . 3037e-01 


-3 . 3 2 50e-01 


-1 . 0638e+00 


-5 . 2 515e-01 


-7 . 9197e-02 


-1. 1663e-01 


1 . 22 55e+00 


-2 . 7079e-01 


2. 6852e+00 


5. 4947e-01 


2 . 160 le-01 


1 . 5287e-01 


8 . 0992e-02 


1 . 87 17e-0 1 


2.354 le-01 


1 . 9563e-01 


1 . 5877e-02 


4 . 0291e-02 


1 . 4 006e-01 


-1. 6166e-01 


2. 6995e-01 


-3 . 87 51e-02 


2 . 3666e-02 


-1 . 7713e-03 


-1 . 9537e-01 


1 . 14 18e-01 


-4 . 0801e-01 


-3 . 2428e-02 


-3 . 3881e-02 


-1. 3917e-02 


-1. 5651e-02 


2 . 682 4e-03 


-3 . 4 512e-02 


-7 . 6200e-03 


-2 . 7647e-03 


-2 . 0665e-03 
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Columns 13 through 18 



-4 . 8 7 87 e-04 


-4 . 7 600e — 04 


-1. 1648e-05 


3 . 4265e+00 


2 . 609 3e-0 1 


7 . 0 03 3e-01 


-1. 1798e-03 


-4 . 6223e-04 


-5 . 063 6e-05 


3 . 2032e-02 


3 . 9424e+00 


4 . 6 16 3e+00 


-1 . 864 le-01 


-1.778 5e-0 1 


-4 . 58 1 5e-03 


1 . 7292e+01 


6 . 8 088e+00 


9 . 9599e+00 


6 . 2 164e-0l 


3 . 9705e-01 


2 . 167 6e-02 


9 . 4 193e+00 


-1. 1360e+01 


-1. 2205e+01 


-1 . 908 6e+00 


-1 . 594 9e+00 


-5.4 329e-02 


-1 . 294 2e + 01 


-3.292 le + 00 


-5.34 3 8e+00 


-1.9215e-01 


9.8820e-02 


-1 . 3922e-02 


5 . 8151e-01 


7 . 0611e+00 


8 . 3284e+00 


-7 . 37 4 2e-0 1 


-2 . 4 886e-01 


-3 . 2950e-02 


8 . 4 62 5e-0 1 


7.77 82e+00 


9 . l980e+00 


1 . 4 035e + 00 


8 . 0952e-01 


5. 1750e-02 


- 5 . 84 24 e-01 


5. 6082e+00 


6 . 4942e+00 


-5 . 7 313e+00 


-1 . 7434e+00 


-2 . 6308e-01 


-7 . 8856e+00 


1 . 7695e+00 


1. 1612e+00 


2 . 3404e+00 


3 . 5305e+00 


1. 1476e-02 


-1 . 51 54e+00 


2 . 0088e-01 


6 . 03 22e-02 


-8 . 154 3e+00 


-2 . 59 54e + 00 


-3 . 6511e-01 


1.393 le-01 


-1.3278e+00 


-1 .537 5e+00 


2.464 7e+00 


1. 5106e+01 


-1. 0758e+00 


-1 . 34 4 2e-01 


7 . 2827e-01 


8 . 3657e-01 


-4 . 563 le+01 


-6 . 34 4 Oe+Ol 


-3 . 3 130e+00 


-1 . 3589e-01 


-2 . 87 18e-02 


-4.92 66e-02 


3 . 3 3 96e+0 1 


-7 . 13 15e+01 


-6 . 2957e+00 


2 . 869 4e-02 


-1 . 7 3 19e-01 


-1.9933e-01 


-2 . 7014e+00 


-3 . 0793e+00 


-1 . 44 18e+02 


-3 . 84 32e-03 


-7 . 298 le-03 


-8 . 9 8 18e-03 


-4 . 607 8e-02 


-3 . 5124e-02 


-1 . 4 208e-03 


1 . 22 50e-01 


2 . 2 3 08e+00 


8 . 2904e-02 


-1 . 8895e-01 


-6 . 5 13 5e-02 


-8 . 3992e-03 


2 . 5142e+00 


-6. 6192e-01 


-1 . 7140e+00 


1 . 8667e-01 


3.40 2 le-02 


9 . 2 87 2e-03 


1. 0453e-01 


1. 5308e+00 


-1. 4431e+01 


-9.4784e-02 


-3 . 0659e-02 


-4 . 279 le-03 


1 . 4 115e-01 


-5 . 66 14e-01 


1 . 202 3e+01 


-2. 1934e-02 


-4 . 5366e-02 


2 . 5787e-04 


-2 . 7934e-01 


2 . 07 13e-01 


-7 . 7357e+00 


4 . 304 2e-02 


6 . 2808e-02 


3 . 4892e-04 


2.74 02e-0 1 


-2 . 0757e-01 


7 . 78 28e+00 


5 . 4 190e-02 


4 . 7 124e-03 


2 . 864 4e-03 


-1 . 2447e-01 


4 . 5133e-01 


-9 . 8684e+00 


4.763 5e — 02 


1. 5922e-02 


2.13 38e-03 


-5 . 093 le-02 


2 . 79l9e-01 


-5 . 4 5 53e + 00 


1.894 3e-02 


2 . 3352e-02 


2 . 9346e-04 


8.37 2 5e-02 


-5 . 0725e-02 


2 . 19 14e+00 


9 . 9162e-03 


-9. l785e-03 


8 . 5161e-04 


-8.43 2 2e-02 


1 . 8 04 2e-0 1 


-4 . 34 99e+00 


-1. 0099e-02 


-1 . 02 60e-02 


-2 . 2788e-04 


-3 . 62 2 6e-02 


5 . 828 2e-03 


-6 . 0937e-01 


7 . 988 5e-03 


5.104 le-03 


2.78 4 6e- 04 


4 . 8 1 50e — 0 3 


2 . 2 795e-02 


-2 . 96l0e-01 


-5. 3051e-03 


-1 . 8796e-03 


-2 . 34 17e-04 


4 . 2518e-03 


-3 . 04 01e-02 


5 . 84 08e-01 


-8. 5947e-05 


1.537 4 e-04 


-9 . 8 005e-06 


2.9774 e-04 


-2 . 2 13 3e-03 


5. 6798e-02 



Columns 19 through 24 



-3.64 2 6e-01 


-2 . 1595e+00 


2 . 4683e-01 


-2 . 2989e-01 


-2 . 2 162e-01 


1 . 8 63 le-03 


1 . 8269e+00 


-1. 1359e+00 


1 . 9499e+00 


-1 . 562 6e + 00 


8.9748e-01 


-6 . 5908e-01 


7 . 1306e-01 


-1 . 24 53e+01 


3 . 9605e+00 


-3 . 3354e+00 


1 . 3 54 5e-01 


-9 . 093 3e-01 


-6. 6120e+00 


-2 . 5 162e+00 


-5 . 29 15e+00 


4 . 1514e+00 


-3 . 3666e+00 


2 . 02 55e + 00 


3 . 0422e-01 


8 . 8093e+00 


-2 . 0725e+00 


1 . 7 8 19e+00 


3 . 104 le-01 


3 . 7883e-01 


3 . 1979e+00 


-2 . 3 53 5e + 00 


3 . 5104e+00 


-2 . 8 18 le+00 


1 . 564 4e + 00 


-1. 1735e+00 


3 . 4936e+00 


-2 . 7 178e+00 


3 . 87 4 0e + 00 


-3 . 1119e+00 


1. 7064e+00 


-1 . 2899e+00 


2 . 688 le+00 


-1. 2326e+00 


2 . 7 52 2e+00 


-2 . 1996e+00 


1 . 3284e+00 


-9 . 4 593e-01 


1 . 93 93e+00 


4 . 2986e+00 


6 . 0 3 52e-01 


-4 . 0962e-01 


1 . 0512e+00 


-4 . 0076e-01 


3 . 08 03e-01 


8 . 654 6e-0l 


4 . 7 18 4e— 02 


-2 . 3598e-02 


1 . 7023e-01 


-5.373 3e-02 


-6 . 3660e-01 


2 . 9124e-01 


-6. 5159e-01 


5. 2077e-01 


-3 . 1460e-01 


2 . 2 398e-01 


3 . 5737e-01 


-1 . 2442e-01 


3 . 553 8e-01 


-2 . 834 8e-01 


1 . 7731e-01 


-1 . 2 361e-01 


5 . 9127e-03 


9 . 084 4e-02 


-1 . 8 869e-02 


1. 6393e-02 


4 . 59 5 4 e-0 3 


2 . 9990e-03 


-8 . 4524e-02 


3 . 158 le-02 


-8 . 4626e-02 


6. 7536e-02 


-4 . 1897e-02 


2 . 9354e-02 


-2 . 84 59e-03 


4 . 4059e-03 


-3 . 7 398e-03 


3 . 0 324 e-03 


-1. 3509e-03 


1. l698e-03 


-3 . 17 3 le-02 


-1 . 7755e-01 


1. 8591e-02 


-1 . 74 28e-02 


-1 . 9 17 6e-02 


8 . 0964e-04 


-6 . 6 508 e-0 1 


2.854 Oe-0 1 


-6.78 1 6e-0 1 


5.4144e-01 


-3. 3059e-01 


2 . 3398e-01 


-1. 1960e+01 


8 . 167 le+00 


-1.227 8e + 01 


9 . 8902e+00 


-5 . 4 14 8e + 00 


4 . 058 le+00 


-3 . 27 52e + 01 


2 . 54 03e + 01 


-6 . 524 4e+01 


4 . 38 19e+01 


-3 . 6204e+01 


2 . 3 589e+01 


1.247 6e + 01 


-3 . 3664e+01 


3 . 1172e+01 


-3 . 4986e+01 


2 . 4 95 le+01 


-1 . 9603e+01 


-9 . 8155e+00 


6 . 677 4e+01 


- 1 . 1499e+02 


2.394 8e+02 


-1.353 le+02 


8.537 le+01 


4 . 1919e+01 


-3 . 7902e+01 


-6 . 7 84 3e+01 


-9 . 1153e+0 1 


1 . 7 52 le+02 


-8 . 2929e+01 


3 . 6204e+01 


-2 . 3385e+01 


9 . 0979e+01 


-1.732 3e + 02 


-7 . 34 40e+01 


1 . 0393e+02 


-8.9347 e + 00 


2 . 7 2 18e+01 


-8.519 5e+01 


5 . 6306e+01 


6 . 6647e+00 


-1 . 4 028e+02 


2 . 2762e+01 


-3 . 0465e+01 


8.6383e+01 


-9 .337 le+01 


-8. 160 le+01 


1. 3847e+02 


1 . 7222e+00 


-9 . 7 024e+00 


2 . 9 393e+01 


-1.4 3 57e + 0 1 


-4 . 7599e+00 


9 . 1385e+01 


2 . 393 6e+00 


1 . 6034e+00 


-6 . 09 56e+00 


-4 . 3450e-f00 


-1. 1412e+01 


-2.37 50e+0 1 


-3.4 64 3e+00 


1 . 9987e+00 


-4 . 4 4 09e + 00 


1.0620e+01 


1. 5786e+01 


-1 . 24 I5e+01 


-2 . 8422e-01 


4 . 057 5e-01 


-1. 1027e+00 


1. 1095e+00 


1. 2718e+00 


-3 . 7835e+00 
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Columns 25 through 29 



7 . 6382e-01 
6 . 4 187e-01 
4 . 7395e+00 
1.5408e-01 
-3 . 2565e+00 
1.2600e+00 
1. 4313e+00 
7 . 804 le-01 
-1 . 37 60e+00 
-2 . 8 68 5e-01 
-1 . 8 4 57e-01 
8 . 9 0 10e-02 
-3 . 3252e-02 
-2. 1857e-02 
-1 . 98 53e-03 
6 . 2 396e-02 
-1 . 8649e-01 
-4 . 2170e+00 
-2 . 1484e+01 
3 . 1664e+01 
-8. 5469e+01 
8 . 7428e+01 
-7 . 3 125e+01 

1 . 4 883e+02 
-2 . 1400e+02 
-8 . 9390e+01 
-8 . 0190e+00 
4. 3488e+01 
5.743 le+00 



-1. 6414e-01 
2 . 3732e-01 
-4 . 954 3e-01 
-1. 1834e+00 
4 . 8012e-01 
3 . 97 53e-01 
4.27 10e-01 
3 . 7 067e-01 
5 . 214 4e-01 
9. 1768e-02 
-8.78 15e-02 
5.12 10e-02 

5 . 3 962e-03 
-1 . 2 007e-02 
-2 . 4 059e-04 
-1. 3781e-02 
-9 . 3 08 4 e-02 
-1.4053e+00 
-8 . 7977e+00 

4 . 0658e+00 
-2 . 7 7 50e + 01 
2 . 7088e+01 
-3 . 9971e+01 

8 . 4 695e+01 
-6 .4387 e+0 1 
-1. 1436e+02 

8.5858e+01 
9 . 82 15e+00 

7 . 04 60e+00 



-1 . 7 64 0e-01 
-5 . 7108e-02 
-9. 6751e-01 
-3 . 149 le-01 
6 . 9870e-01 
-1 . 2870e-01 
-1. 5214e-01 
-4 . 949 le-02 
3 . 7259e-01 
7 . 3 560e-02 
1 . 1669e-02 
-3 . 47 55e-03 
7 . 2 542e-03 
9.914 5e-04 
2.964 6e-04 
-1.4500e-02 
1. 0728e-02 
4 . 13 95e-01 
1. 5914e+00 
-4 . 5088e+00 
7 . 44 58e+00 
-7 . 37 69e + 00 
7 . 1874e+00 
-2 . 14 2 3e+0 1 
5.727 4 e + 0 1 
-3 . 152 3e+01 
—4 . 04 64 e+0 1 
3 . 192 5e+01 
-2 . 93 37e-01 



1. 7079e-01 
2.472 le-02 
8.94 17 e-0 1 
3 . 9862e-01 
-6 . 58 6 le-0 1 
7 . 0 1 69e-02 
8.74 56e-02 
4 . 0556e-03 
-3 . 7915e-01 
-7 . 3678e-02 
-9. 1201e-04 
-2 . 3655e-03 
-6.88 13e-03 
4 . 009 3 e-04 
-2 . 3268e-04 
1.407 0e-02 
4 . 611 6e-04 
-2. 1229e-01 
-4 . 6398e-01 
3.4961e+00 
-3. 5528e+00 
3 . 6175e+00 
-1. 8422e+00 
1 . 0333e+01 
-4 . 14 14e+01 
9. 1100e+00 
6. 0359e+01 
-1. 3722e+02 
-2 . 1276e+01 



bcgf = 

5 . 2502e+00 
-3 . 2830e-01 
2 . 597 0e + 0 1 
1.5589e+01 
-1 . 9609e+01 
2 . 1893e-01 
5 . 578 3e-01 
-1.4366e+00 
-1 . 2 3 10e+01 
-2 . 3 52 2e + 00 
3 . 4 165e-01 
-2 . 7 67 Oe-0 1 
-2 . 064 6e-01 
6 . 0765e-02 
-5 . 2l77e-03 
-5 . 7008e-02 
3 . 0829e-0 1 
-5 . 129 6e-0 1 
1. 6850e-01 
-2. 2395e-01 
2 . 593 5e-01 
-1.84 38e-01 
-8. 1148e-02 

8 . 4 666e-02 
-1 . 0273e-01 
-3. 0092e-02 

3 . 1136e-03 
8. 3068e-03 

1 . 4 001e-03 



6 . 4 618e-01 
2 . 9101e+00 
7 . 3091e+00 
-7 . 1255e+00 
-4 . 1 4 07e+00 
5 . 2816e+00 
5 . 84 52e + 00 
4 . 0563e+00 
2 . 5979e-01 
-5 . 2 653e-02 
-9 . 6027e-01 
5. 1899e-01 
-3 . 9 167e-02 
-1. 2386e-01 
-5 . 8883e-03 
-1 . 3 306e-02 
7. 1960e-02 
-1. 1973e-01 
3.933 le-02 
-5.227 4e-02 
6 . 053 6e-02 
-4 . 3038e-02 
-1 . 894 le-02 
1. 9762e-02 
-2 . 3978e-02 
-7 . 02 39e-03 
7 . 2 677e-04 
1 . 9389e-03 
3 . 2679e-04 



7 . 1694e-03 

9 . 4 652e-03 
2 . 2894e-02 
4 . 8928e-02 
2 . 1497e-02 
1 . 5760e-02 
1 . 6892e-02 

1 . 4 898e-02 
2 . 2 23 3e-02 
3 . 93 59e-03 
3 . 5297e-03 
2 . 0680e-03 
2 . 3989e-04 
4 . 8 4 34e-04 
8.907 5e-06 
6. 0100e-04 

3 . 7 4 62e-03 
5.584 le-02 
3 . 5200e-01 
1 . 4919e-01 
1. 0843e+00 
1 . 0655e+00 
1. 6118e+00 
4 . 0479e+00 
3 . 2024e+00 
9 . 2255e+00 
9 . 2992e-0 1 

4 . 07 8 3e-»-00 
1 . 3857e + 02 
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ccgf = 

Columns 1 through 6 

0 

0 

Columns 7 through 12 



Columns 13 through 18 



0 -6 . 4236e-01 1.1544e-01 

0 -8.347 4 e-02 -1.3417e+00 



6 . 0999e-02 
-1. 5794e+00 



Columns 19 through 24 



1 . 44 64e-01 
-6. 1148e-01 



3 . 582 9e-01 
4 . 3076e-01 



3.497 0e-02 
-6.6611e-01 



-2 . 19 9 6e-02 7 . 9 072e-0 2 

5.3449e-01 -2 . 9948e-01 



-2.7843e-02 
2 . 2334e-01 



Columns 25 through 29 

-1. 1672e-01 4 . 07 3 6e-02 

-2 . 337 3e-01 -7.69 55e-02 



3 . 07 55e-02 
2 . 3094e-02 



-3 . 1054e-02 
-1. 19 98e-02 



-1 . 74 17e-03 
3 . 0561e-03 



dcgf = 
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APPENDIX E 



Poles of The 3 I/O Reduced Order Closed Loop System 



-9. 1468e+03 
-5 . 5159e+02 + 
-5 . 5159e+02 - 
-2 . 2777e+02 + 
-2 . 2777e+02 - 
-1 . 4457e+02 
-1 . 4319e+02 
-1 . 1949e+02 
-8 . 9079e+01 + 
-8 . 9079e+01 - 
-1. 7935e+01 + 
-1 . 7935e+01 - 
-5 . 2531e+01 + 
-5 . 2531e+01 - 
-4 . 9021e+01 + 
-4 . 9021e+01 - 
—3 . 1654e+01 + 
-3 . 1654e+01 - 
-2 . 9006e+01 
-2 . 0127e+01 
-8 . 7118e+00 + 
-8 . 7118e+00 - 
-8 . 094 5e+00 + 
—8 . 094 5e+00 - 
-1. 4512e+01 
-1 . 8891e+00 
-1 . 9503e+00 
-1 . 0702e+01 
-2 . 7 168e+00 
-8 . 6993e+00 



5 . 3384e+02i 
5 . 3384e+02i 
2 . 3189e+02i 
2 . 3189e+02i 



1. 8937e+01i 
1 . 8937e+01i 
6. 0318e+01i 
6. 0318e+01i 
4 . 8268e+01i 
4 . 8268e+01i 
4 . 9207e+0li 
4 . 9207e+01i 
9 . 6963e+00i 
9 . 6963e+00i 



9 . 5608e+00i 
9 . 5608e+00i 
6. 1691e+00i 
6. 1691e+00i 



75 



REFERENCES 



1. Doyle, J. C. and Stein, G., Multivariable Feedback Design: Concepts for a 
Classical/Modern Synthesis, IEEE Transactions on Automatic Control, Vol. AC- 
26, No. 1, pp. 4-16, February 1981. 

2. Safonov, M. G., Laub, A. J., Hartman, G. L., Feedback Properties of 
Multivariable Systems: The Role and Use of The Return Difference Matrix . 
IEEE Transactions on Automatic Control, Vol. AC-26, No. 1, pp. 47-65, February 
1981. 

3. Lehtomaki, N. A., Practical Robustness Measures in Multivariable Control 
System Analysis, Ph.D. Dissertation, Massachusetts Institute of Technology, 
Cambridge, Massachusetts, May 1981. 

4. Gordon, N. C., Utilization of Numerical Optimization Techniques in the 
Design of Robust Multi-Input Multi-Output Control Systems . PhD. 
Dissertation, Naval Postgraduate School, Monterey, California, September 1984. 

5. Chiang, R. Y., and Safonov, M. G., Robust-Control Toolbox. User’s Guide, pp. 
R4-R33, RR9, RR10-RR11, RR26-RR30, The Mathworks, Inc., June 1988. 

6. Postlethwaite, I., and others, H„ Control System Design: A Critical 

Assessment International Federation on Automatic Control, Vol. 8, pp. 328- 
333, 1987. 

7. Kwaakemaak, H., and Sivan, R., Linear optimal Control Systems, John Wiley 
and Sons, Inc., 1972. 

8. Ogata, K., Modern Control Engineering. Prentice-Hall, Inc., 1970. 

9. Safonov, M. G., and Chiang, R. Y., CACSD Using the State Space L_ Theory: 
A Design Example . IEEE Transactions on Automatic Control, Vol. 33, No. 5, 
pp. 477-479, May 1988. 

10. Jane’s All the World’s Aircraft 1987-1988 . Jane’s Information Group Limited, 
1988. 

11. V. C. Gordon and D. J. Collins, Multi-Input Multi-Output Automatic Design 
Synthesis for Performance and Robustness . AIAA Journal of Guidance, 
Control, and Dynamics, Vol. 9, No. 3, pp. 281-287, May-Jane 1986. 

12. Hidenori Kimura, Robust Stabilizabilitv for A Class of Transfer Functions . 
IEEE Trans. Automat. Contr., VOL. AC-29, No. 9. pp. 788-793. Sept. 1984. 



76 



13. W. L. Rogers, Applications of Modern Control Theory Synthesis to A Super- 
Augmented Aircraft . M.S. Thesis, Naval Postgraduate School, Monterey, 
California, June 1989. 

14. Bernard Friedland, Control System Design: An Introduction to State Space 
Methods . McGraw-Hill Book Company, 1986. 

15. Bruce A. Francis and others, H_ - Optimal Feedback Controllers for Linear 
Multivariable Systems . IEEE Trans. Automat. Contr., vol. AC-29, No. 10, pp. 
888-899, Oct. 1984. 



77 



INITIAL DISTRIBUTION LIST 



No. copies 

1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22304-6145 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, California 93943-5000 

3. Chairman, Code 69Hy 2 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93943-5000 

4. Professor D. J. Collins, Code 67Co 5 

Naval Postgraduate School 

Monterey, California 93943-5000 

5. Professor Liang-Wey Chang, Code 69Ck 1 

Naval Postgraduate School 

Monterey, California 93943-5000 

6. Main Library 1 

P. O. Box 8243, Ta-Hsi, Tao-Yuan, 33500 

Taiwan, Rep. of China 

7. Mr. Ta-Chieh Hsu 3 

P. O. Box 8243-23, Ta-Hsi, Tao-Yuan, 33500 

Taiwan, Rep. of China 

8. Mr. Tor W. Jensen Code 6013 1 

Naval Air Development Center 

Warminster, Pennsylvania 18974 

9. Mr. Joe Gera 1 

NASA Dryden Flight Research Center 

P.O. Box 273 
Mail Code OFDC 
Edwards, California 93523 



78 



10. Mr. Thomas Momiyama 1 

AIR 931 

Naval Air Systems Command 
Washington, D.C. 20361-0001 

11. Mr. George Derderian 1 

AIR 931E 

Naval Air Systems Command 
Washington, D.C. 20361-0001 



79 



Ic/C 



Thesis 
H82973 Hsu 

cA Application of H ot) 

method to modern fighter 
configuration. 



Thesis 
H82973 Hsu 

Application of H 00 
method to modern fighter 
configuration. 



